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Chapter 1 
Introduction 



1.1 Reasons for gauge theory of the strong interactions 

Quantum Chromodynamics (QCD) is an SUc{S) colour gauge theory, describing interac- 
tions of quarks in the fundamental representation. We begin by summarizing the standard 
arguments in favour of QCD as the theory for strong interaction physics (see for instance 



1. At an empirical level, light hadrons (mesons and baryons) can be classified according 
to representations of SUf{3) (the subscript / indicates flavour). It is observed that 
precisely those representations of SUf{3) are physically realized that may be obtained 
by decomposing the direct products 3 3* (mesons) or 3 ® 3 (S) 3 (baryons) as direct 
sums of irreducible representations of SUf{3). 



It is an important fact that the fundamental representations of SUf{3), 3 and 3*, are 
not among the physically realized representations. 

2. This leads to the hypothesis of quarks and antiquarks as the fundamental constituents 
of hadrons (Quark Model of hadrons). One finds that there are three "fiavours" of light 
quarks (antiquarks), which transform according to the fundamental representation 
3 (3*) of SUf{3). These three flavours are referred to as "up"(m), "down"((i) and 
"strange"(s), and they have electric charges (f ? ~|) ^|) respectively. Consideration 
of SUf{3) as an approximate symmetry of the strong interactions leads to the above 
mentioned classification of hadrons. The approximate nature of SUf{3) is due to 
the fact that the three quark fiavours have different masses and electric charges. 
In the Non-relativistic Quark Model, the intrinsic angular momenta associated with 
hadrons are obtained by combining S'f/spin(2) with SUf{3). This leads to a reasonably 
successful SU{6) classification scheme. The orbital excitations are then included by 
classifying hadrons under SU{6) ® Ospace(3). 
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3. Absence of quarks among the observed hadrons requires them to exist only as con- 
stituents of bound states. The quarks are assigned colour charge, and the strong 
interactions bind them into colour neutral hadrons. This is the phenomenon of con- 
finement. In the Non-relativistic Quark Model, a meson is a bound state of a quark 
and an antiquark, while a baryon is a bound state of three quarks. 

4. Mesons are bosons, while baryons are fermions. This requires quarks to be fermions, 
and baryon wavefunctions to be fully antisymmetric when all their degrees of free- 
dom are taken into account. In the Non-relativistic Quark Model, the lowest lying 
baryons are in zero relative angular momentum states, and their wavefunctions are 
symmetric in spin-flavour degrees of freedom. For example, the wavefunction of the 

A++ (l"*") I'"T'"T'"t)i where the up-arrow (j) denotes the spin Sz = for each 
quark. These wavefunctions are reconciled with Fermi statistics by postulating a 
new internal quantum number for quarks, i.e. colour. If each quark has three, other- 
wise indistinguishable, colour states, which transform according to the fundamental 
representation oi SUc{S) (subscript c stands for colour), then Fermi statistics is saved 
by using a totally antisymmetric colour wavefunction for baryons. 

5. This assignment makes strong interactions invariant under global SUc{3) transfor- 
mations. The states may then be classified according to various irreducible repre- 
sentations of SUc{3). Using the fact that physical hadrons are colour neutral, i.e. 
singlets under SUc{3), we can understand why only qq and qqq states (and not q, qq, 
qqqq etc.) exist in nature — the singlet representation appears only in the 3® 3* and 
3 (8) 3 (8) 3 products. 

6. Data from deep inelastic e—p scattering (which involves a large momentum transfer 
to the proton via the exchange of a virtual photon) suggest that the proton can be 
understood as a bound state of constituents called "partons", which interact very 
weakly with each other at high energies. At low energies, the interactions among 
the partons becomes strong and they are strongly bound to each other. These ideas 
form the basis of Feynman's Parton Model. Non-abelian gauge theories can provide a 
field-theoretic understanding for Parton Model, because it has been shown that only 
non-abelian gauge theories are asymptotically free, i.e. their interactions vanish as 
energy E ^ oo but become strong at a scale E ~ Aqcd (Aqcd denotes the O(GeV) 
confinement scale). 

7. Experimentally, there is no evidence for any flavour dependence of strong interactions; 
it is found that flavour-dependent effects can be explained by quark mass and electric 
charge difi'erences. We have already argued that strong interactions should have an 
exact global SUc{3) symmetry. Converting this to a local symmetry, SUc{3) becomes 
the non-abelian gauge theory of strong interactions. 

8. The concept of three colour states for quarks is supported by at least two other 
powerful experimental results. 
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(a) The parton model cross-section calculation for the process e'^ + e — > hadrons: 
At large energies {E ^ Aqcd), the ratio 

^ a {e^e~ — > hadrons) 

is predicted to be 

R^ZY^QI = 3x 

= — (including quarks up to h) 

The experimental value of R is in good agreement with this prediction and in 
poor agreement with the colourless prediction y. 

(b) The decay rate for 7r° 27: This is again proportional to the number of colours, 
and the observed decay rate agrees with three colour states for the quarks. 

All this analysis suggests that the theory of strong interactions is a non-abelian gauge the- 
ory based on the gauge group SUc{'i) with quarks in the fundamental representation. This 
theory is precisely defined with a single gauge coupling and quark masses as its parame- 
ters. For its confirmation, it must quantitatively explain quark confinement, spontaneous 
breakdown of chiral symmetry (i.e. light pions), the hadron spectrum, partial decay widths, 
cross-sections, form-factors, structure functions, and so on. 
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1.2 QCD 

QCD is a non-abelian gauge theory based on the gauge group SUc{'i). We assume that 
there are Nf quark flavours (experimentally Nf = 6). Each quark of a specific flavour n 
comes in 3 colours, and transforms according to the fundamental representation oi SUc{^). 
The QCD Lagrangian for Nf flavours of quarks is given by 

n=l 

Various symmetries of the QCD follow from C. 

• Invariance imder the (3 + l)-dim Poincare group, which consists of time translations, 
space translations, spatial rotations, Lorentz boosts, and their combinations. Con- 
servation of energy, linear momentum and angular momentum follow as a result of 
these space-time symmetries. 

• Charge conjugation, Parity and Time-reversal invariance. 



CHAPTER 1. INTRODUCTION 



SUc{S) colour gauge invariance. 

Invariance under ^„ — > exp [iOn] and ^„ — > exp [— i^n] ■ This corresponds to 



U{1) U{1) ® • • • ® U{1) = [U{l)f^' 



Nf times 



global invariance, which implies conservation of the number of quarks minus the 
number of antiquarks for every flavour. 

• When all = 0, invariance under 

*n ^ (exp iOktk + iOktk75\) *m , *n ^ *m (exp -iOktk + iOktk^bl) 

\ L _ _ J / nm \ L _ _ J / mn 



exp 



ie + i^75j , ^ *n exp ^-ie + ^^75] 



where tfc are the Nj — 1 matrices generating the Lie algebra of SU (Nf) group, and 

9k and 6k are arbitrary real parameters (see Appendix A). Thus the global [f/(l)]^-^ 
gets extended to SUl (Nf) SUr {Nj) Uvil) ® which is known as the the 

chiral group of symmetries. The Uv{l) symmetry is exact and corresponds to the 
conservation of Baryon number. The C/a(1) symmetry is anomalous, i.e. although 
it is a symmetry of the Lagrangian, it is not a symmetry of the quantum dynamics. 
Since parity doubling is absent in the observed hadron spectrum, one expects the 
SUl {Nf) ^SUr {Nf) symmetry to be spontaneously broken to its diagonal subgroup 
SUv{Nf). A signature of this spontaneous symmetry breaking is provided by the 
non-zero expectation value of the chiral condensate, 7^ 0, and Nj — 1 massless 

Goldstone bosons. (For Nf — 2, the triplet of pious are identified as the pseudo- 
Goldstone bosons). 

The axial (^4) part of the chiral symmetry is broken explicitly by the non-zero bare 
quark masses, while the vector {V) part is explicitly broken due to the fact that the 
bare quark masses are unequal. Thus the mass terms account for explicit flavour 
symmetry breaking effects in QCD. 



An important aspect of 4-dim QCD which arises as a result of the non-abelian nature of 
the gauge group SUc{3) is asymptotic freedom. The coupling constant goes to zero as the 
energy scale at which it is defined is increased. Thus the behaviour of the Green's functions 
when all the momenta are scaled up by a common factor is governed by a theory where 
g ^ 0. For processes which involve large momentum transfers, perturbation theory in 
powers of g can be reliably used, and it is in this domain that QCD has been extensively 
tested. Factorization of long distance physics (which is process independent) from the 
short distance distance physics (which is process specific) gives perturbative QCD a clear 
predictive power. Perturbation theory also indicates that ^ — > 00 at a low energy scale, 
typically labeled Aqcd- The strong coupling region {g very large) has been studied for 
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QCD formulated on a space-time lattice using high temperature expansions. Various non- 
perturbative features of hadronic physics (e.g. confinement, chiral symmetry breaking and 
quantum numbers of light hadrons) emerge naturally when g is large. To connect the two 
regions, QCD has to be analysed at intermediate values of g. It has not been possible to do 
this analytically so far. Numerical simulations have extensively explored this intervening 
region, however, and they show that the region of large gauge coupling is analytically 
connected to the region of small gauge coupling, without any phase transition in between. 
This implies that the symmetries of the theory are the same in both the strong coupling 
and the weak coupling regions. 

1.3 Summary of the thesis 

There is little doubt that QCD is the theory behind hadronic physics. Our main difficulty 
lies in performing quantitative calculations using it. Such calculations are important be- 
cause they test the success of QCD in low energy domain. For example, the cross-section 
for a process involving hadrons, in general, has a sizable non-perturbative component re- 
lated to hadronic substructure. If such processes are to be used to identify physics beyond 
the Standard Model, then it is important to have a good theoretical control over the non- 
perturbative contribution. 

The study of bound states in QCD is made difficult due to [221 

• the unlimited growth of the running coupling g in the infrared region, which invali- 
dates use of perturbation theory; 

• confinement, which requires potentials that diverge at long distances as opposed to 
the Coulomb/ Yukawa potentials used in perturbation theory; 

• spontaneous chiral symmetry breaking, which cannot be demonstrated using pertur- 
bation theory; and 

• the non-perturbative structure of the QCD vacuum. 

Before one can quantitatively understand QCD, it is necessary to establish a suitable 
framework for studying these non-perturbative effects. A lattice formulation is a natural 
way to define QCD so that all the non-perturbative physics is accessible. Much progress 
has been made since QCD was formulated on a space-time lattice by Wilson in 1974. 
In this thesis, we study QCD formulated on a transverse lattice, which combines techniques 
from light-front field theory and lattice gauge theory. An approximation scheme based on 
l/Nc expansion and transverse strong gauge coupling expansion allows us to systematically 
study this transverse lattice QCD. 

The use of light-front framework in studying bound states has several advantages over 
the equal-time framework. In this framework, the theory is evolved along a light-like 
direction (x^). The light-front wavefunctions, which encode the bound state structure, are 
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the eigenstates of the light-front Hamiltonian P~ . Since the longitudinal and transverse 
boosts are purely kinematical, it becomes easy to separate the centre-of-mass motion from 
the internal motion. Furthermore one can choose to work with internal variables which 
are manifestly boost invariant. The computation of the light-front wavefunctions leads 
to estimates of a wide range of observables. Implementing a cut-off that removes degrees 
of freedom with exactly zero longitudinal momentum, one can force the vacuum to be 
trivial. The bound state wavefunctions can then be expanded in the light-front Fock space 
of its constituents. In the large-A^'c limit, only the lowest Fock space amplitudes survive. 
The non-dynamical longitudinal gauge field provides a linear confining potential in the 
^"-direction (see Appendix B). 

The lattice cut-off in the transverse directions provides a non-perturbative gauge- 
invariant regularization, albeit at the cost of rotational invariance. In the strong transverse 
gauge coupling limit, there is a linear confining potential in the transverse directions. We 
expect QCD formulated with two continuum (light-front) and two lattice (transverse) di- 
mensions to be closer to real continuum QCD, than the extensively studied version with 
all the four dimensions on lattice. In the limit of \arge-Nc and strong transverse gauge cou- 
pling, the dynamics of QCD along the light-front and the transverse directions essentially 
factorizes, and the theory can be solved in a closed form. This approach to studying QCD 
was outlined by Patel j22|. Here we present the explicit solution, and that is the original 
contribution of this thesis. 

In Chapter 2, we formulate QCD in the large- A^'c and strong transverse coupling limits, 
and then exactly integrate out all the gauge degrees of freedom to obtain the generating 
functional for quark-antiquark bilinears. We study the chiral properties of our limiting 
theory in Chapter 3, for naive as well as Wilson fermions, and obtain a recursive relation for 
the chiral condensate. In Chapter 4, we obtain the homogeneous integral equation satisfied 
by the meson states of our theory. Comparison of this equation with the corresponding one 
for the 't Hooft model allows us to infer many physical meson properties. Our results are 
consistent with phenomenological expectations, and this is the first time that such results 
have been obtained from (3 + l)-dim QCD, with only quark and gluon degrees of freedom 
and no ad hoc model assumptions. Extraction of precise values requires numerical solution 
of the integral equation; we have not yet carried that out and we describe our outlook for 
further investigations at the end of Chapter 4. 

Three appendices supplement our analysis. Our notation and conventions are listed in 
Appendix A. The known results for mesons in the 't Hooft model and in strong coupling 
lattice QCD are rederived, in Appendix B and in Appendix C respectively, using the same 
methodology as followed in the thesis. That allows convenient comparison, as well as 
demonstrates the advantage of our approach over these two well-studied approximations 
to QCD. 
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QCD at large Nq and strong transverse 
gauge coupling 

2.1 Large- A^c limit 

QCD is a non-abelian gauge theory based on the gauge group SU{3). One would like to see 
all the qualitative and quantitative features of hadronic physics (confinement, chiral sym- 
metry breaking, the hadron spectrum, form factors, partial decay widths, cross-sections, 
etc.) emerge from QCD. For most of these properties, the low momentum structure of 
Green's functions is important. But QCD is an asymptotically free theory, with the running 
gauge coupling increasing as the scale at which it is defined decreases. The low momentum 
behaviour of Green's functions, therefore, is governed by a theory with a large coupling, 
making ordinary perturbation theory a completely ineffective tool of computation. 

In such a situation, one has to look for alternative approximation schemes. A promising 
approach is the 1/Nc expansion of a gauge theory based on the SU (Nc) group. The 
subject of the large-iVc limit of SU (iVc) gauge theories with quarks in the fundamental 
representation was initiated by 't Hooft |^, wherein he showed that l/N^ can be treated 
as a small expansion parameter in these theories. The usual perturbative expansion reveals 
very little information about the spectrum of the theory, while the lavge-Nc limit retains 
many important non-perturbative aspects of the theory. 

Witten showed that many qualitative features of hadronic physics can be under- 
stood in the framework of the l/Nc expansion. In particular, if one assumes confinement, 
the large- A'^c limit produces a weakly interacting theory of an infinite number of mesons 
with coupling of the order of l/N^. Baryons are solitons of the theory, with mass propor- 
tional to Nc, and with size and shape A^^c-independent. Let us now have a quick look at 
some aspects of the large- A'c limit of QCD. This section is largely based on reviews by Cole- 
man Witten |3I3, Manohar [TH| and Makeenko [12J, and follows the power-counting 
conventions of Witten. 

A sensible large- A'c limit of any quantity in QCD is obtained as follows: 

• Generalize the colour gauge group from SU{3) to SU (Nc), 
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• holding g fixed, take large- A^^ limit of each Feynman diagram contributing to that 
quantity, and sum the leading contributions to all orders in g. 

The central idea of the large- A^c limit of QCD is that gluons have many more colour states 
than quarks {N^ — 1 compared to Nc), and so quark radiative corrections are suppressed 
relative to gluonic radiative corrections when one holds g fixed. Also planar type of glu- 
onic radiative corrections are more important (larger powers of Nc) than non-planar type. 
Thus the gluons effectively provide a mean field environment through which the quarks 
propagate. 

In more concrete terms, the Lagrangian density for SU (Nc) gauge field theory with 
quarks in the fundamental representation is given by 

^ = -|^r^;. + ^^ i^^D, - m],^ , 

with 

{D^^),^d^^,-i^^{ta),,Ayj, 

Fr = d^A: - + -^UcA^^A'^c , 

and 

[ta, h]- — ifabctc , tr {tatb) — Sab ■ 

Repeated colour indices are summed here; see Appendix A for notations and conventions. 
We have left out gauge-fixing terms and renormalization counter-terms. Inclusion of these 
terms is trivial and does not alter the nature of qualitative results obtained in their absence. 
Holding g fixed and taking Nc ^ oo, one can show the following: 

• The leading diagrams in the large- A^c limit are planar diagrams with quarks (if any) 
on boundaries. 

— If Gi represent local gauge invariant operators (which cannot be split further 
into smaller gauge invariant parts) made up exclusively of gauge fields, then 
the n-point Green's function, the connected n-point Green's function and the 
n-point IPI vertex function are of order N^"', N^ and Nc^^~^^ respectively, 

{G,G2...Gr.) oc ATf , 

rgj(i,2,...,n) oc AT^a-). 

— If Bi represent gauge invariant (local or non-local) quark bilinear operators 
(which cannot be split further into smaller gauge invariant parts), then the m- 
point Green's function, the connected m-point Green's function and the m-point 
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IPI vertex function are of order N^, Nc and A*"^ respectively, 

{BiB2...B^) oc NT, 

. . . CX AT,, 

r|'2|(l,2,...,m) oc N]-^. 

— In general, with n-gluonic and m-quark bilinear operators (m > 0), we have 
for the (n + m)-point Green's function, the connected (n + m)-point Green's 
function and the (n + m)-point IPI vertex function, 

{Gl...GnB,... Err.) conn « , 

r|:^;^\(l,...,n;l,...,m) oc nI~'-- . 

• The diagrams with minimum number of quark loops dominate; each internal quark 
loop is suppressed by a factor oil/Nc. 

• Expansion in powers of l/N^ becomes a topological expansion for the diagrams. 
Consider a general diagram contributing to the connected correlation function with 
n-gluonic and m-quark bilinear operators. Any such diagram can be associated with 
a 2-dim surface by associating colour index loops with faces of polygons, gluon prop- 
agators with edges where polygons touch each other and vertices with points where 
corners of different polygons fuse. The contribution of such a diagram is 0{N^), 
where x is the topologically invariant Euler characteristic of the associated 2-dim 
surface. If the 2-dim surface has H handles and L boundaries, then % = 2 — 2H — L. 
[L is the sum of the number of internal virtual quark loops and the number of exter- 
nal boundaries in the diagram.) For a connected {n + m)-point correlation function, 
a maximum of m external boundaries are allowed. The leading contribution arises 
from "planar" diagrams which have genus H = 0, and contribution of diagrams with 
larger genus H or more boundaries L is suppressed. 

• Distinction between gauge groups SU (Nc) and U (Nc) is immaterial at the leading 
order; they differ from each other by O (l/N^) terms. 

It follows that in the large- A'c limit, all gauge invariant correlation factorize; the fluctuations 
are suppressed by powers ofl/Nc and vanish when Nc = oo. This essentially means that the 
functional integral for the theory gets saturated by contribution from a single gauge orbit, 
known as the master orbit and represented by a set of gauge equivalent vector potentials. 
If one knows the master gauge field, then, in principle, one can compute the leading Nc 
contribution to any quantity. In this sense, the large- A'c limit is like a classical limit of the 
theory. 

Assuming that large- A'c QCD is a confining theory, the following qualitative picture of 
glueballs and mesons emerges PO] : 
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• At A^'c = oo, glueballs (mesons) are free, stable and non-interacting. Glueball (meson) 
masses have smooth finite Hmits as Nc ^ oo, and the number of glueball (meson) 
states is infinite. Specifically, if G and B are gauge invariant field operators for a 
meson and a glueball respectively, with appropriate quantum numbers, then 



• With r[;^^j (1,2, . . .n) oc Nc^^ and (0 | G | gluebaU) oc N^, it follows that the am- 
plitude for a process involving n glueballs is proportional to iV^~". In particular, the 
amplitude for a glueball decaying in to two other glueballs is O (1/iVc), and glueball- 
glueball scattering amplitude is 0{1/N^). 

• With rj'^i (1, 2, ... m) oc Nl~"^ and (0 | 5 | meson) oc ^/Wc, it follows that the ampli- 

tude for a process involving m mesons is proportional to Nc .In particular, the 
amplitude for a meson decaying in to two mesons is O (l/i/iVc), and meson-meson 
scattering amplitude is O (l/Nc). 

• More generally, the amplitude for a process involving n glueballs and m mesons 
(m > 0) is of order O (^Nc'""'""^^^ . In particular, the amplitude for a glueball to mix 
with a meson is O [l/^/N^), the amplitude for a glueball to decay in to two mesons 



is O (1/iVc), the amplitude for a meson to decay in to two glueballs is O [l/Nc 



and the glueball-meson scattering amplitude is O (l/N^). 

• Arbitrary n-point Green's functions and n-point scattering amplitudes are given by 
sums of tree diagrams with effective local vertices and propagators corresponding to 
exchange of physical hadrons (glueballs and mesons) . 

This qualitative picture agrees very well with observed phenomenology: 

1. Suppression of the quark sea in hadronic physics and the fact that mesons are ap- 
proximately pure qq states: In large- A^c limit quark loops are suppressed by factors 
of 1/Nc, which is a consequence of the fact that there are O (A^) gluon degrees of 
freedom, but only O (Nc) quark degrees of freedom. Thus, the quark sea is absent at 
A"c = oo- 

2. The absence, or at least suppression, ofqqqq exotics: Mesons are non-interacting at 
Nc = oo, which means that two mesons would not be bound together in to an exotic 
state. Presence of exotics would show up as simple poles in the connected Green's 
function {E{x)E{y))^^^^^, where E{x) = qlqiqjqj^x). But to leading order in 1/Ac, 
{E{x)E{y))^^^^ factorizes to {qiqiix)qkqk{y))l^^^, and we have freely propagating non- 
interacting mesons instead of exotics. 



( I G I gluebaU ) oc N, 
( I I meson ) oc v 
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3. Resonant two-particle final state domination (whenever kinematically allowed) of 
multi-particle decays of unstable mesons: It is observed that meson decays proceed 
mainly through resonant two-particle states. For example, 5(1237 MeV) decays to 
4:71, but mainly via B ujtx, with subsequent oj 37r. (It is generally believed 
that the tendency for decays to be two-particle dominated persists, even when the 
larger phase space of multi-particle decays is accounted for.) The explanation for 
this feature from the point of view of large- iVc limit is that the decay B — > cuvr is 

O (l/v'iVc), while the direct, non-resonant decay B — > txtiiiv: is 

4. Zweig's rule: Constituent quark-antiquark production is suppressed compared to re- 
arrangement of existing constituents inside hadrons. For example, the decay 0(ss) 
K~^{su) + K~{us) has considerably less phase space compared to 4>{ss) — > 'n"^{du) + 
Ti'^iuu — dd) + n'iud). Nevertheless, K^K~ is the dominant decay mode. In 
large- A'^c limit, creation/annihilation of every quark-antiquark pair gives a suppres- 
sion factor oil/Nc. In 0-decay, there is an extra 1/y/Nl suppression factor due to the 
fact that 7r+7r°7r^ is a three-particle decay while K^K~ is a two-particle 
decay. Overall, the tt+tt^tt" decay mode is suppressed relative to the K^K~ decay 
mode by I/nI^"^ . 

5. Meson multiplet structure as nonets of flavour SU {?>)■. If m, d and s quark masses 
were equal, then there would be singlet-octet degeneracy for mesons to leading order 
in 1/Nc. This is because the diagrams that split singlets from octets involve qq 
annihilation, and are of order 

If one restricts QCD to two dimensions (see Appendix B), it is possible to study various 
dynamical features of the theory analytically (many of these features continue to hold good 
in four dimensions). But it is only the Nc = oo limit which has been solved in a closed form, 
't Hooft obtained an integral equation for the meson wavefunction in 2-dim N^ = oo 
QCD, with eigenstates corresponding to an infinite number of bound states. Subsequently, 
Callan, Coote and Gross [101, Einhorn [12], Brower et al. Cardy (H], Einhorn et al. [T3] . 
and Shei and Tsao |24| further explored questions relating to quark confinement, hadron 
states, scattering amplitudes and high energy behaviour. The key features of this solution 
are: 

• Choosing the light-front gauge A"*" = + = eliminates the non-linear self- 
coupling of the gauge fields (gluon trilinear and quartic terms). This gauge is also 
free of ghosts. 

• The dynamics of gauge fields (terms involving A~) can be replaced by an instanta- 
neous linearly confining potential between quark bilinears. 

• The confining potential is defined using a principal value prescription, which keeps 
its Fourier transform infrared finite. 
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• To study mesons in the light-front gauge, one needs to sum only the planar "rainbow" 
and "ladder" diagrams. 

• The choice of light-front coordinates together with the light-front gauge simplifies the 
Lorentz index structure of the theory, and eliminates the non-dynamical components 
in terms of the dynamical ones. 

• The mass term in the quark propagator gets additively renormalized, while the gluon 
propagator and the quark-gluon vertex receive no corrections. 

• The quark-antiquark scattering amplitude shows that there are no continuum states 
in the spectrum; only discrete meson bound states exist. 

• The meson spectrum (masses and wavefunctions) is determined by the eigenvalues 
of an integral equation that characterizes the singularities of the quark-antiquark 
scattering amplitude. It lies on an approximately linear trajectory in the — n 
plane, where n is the radial excitation quantum number. 

• Spin does not exist in two dimensions, but parity is a good symmetry. The lightest 
meson has negative parity, and the subsequent meson states alternate in parity as a 
function of n. 

• In the large-A^c limit, there is a non-zero chiral condensate, ^ 0, even in two 
dimensions. 

Many of these features are explicitly described in Appendix B. Thus (1 + l)-dim large- iVc 
QCD displays confinement and chiral symmetry breaking as anticipated in (3 -|- l)-dim 
QCD, even though the origins of these features are quite different in the two theories. The 
ability to solve (1 + l)-dim large- iVc QCD proves extremely useful, when attempting to 
solve (3 -|- l)-dim QCD on a transverse lattice. 

2.2 Strong coupling limit 

QCD is an asymptotically free theory, so its running gauge coupling grows as its cut-off 
is lowered. Perturbative renormalization group evolution shows that the gauge coupling 
depends logarithmically on the cut-off, and diverges at a scale Aqcd = O(GcV). Strong 
coupling expansion, i.e. an expansion in inverse powers of the gauge coupling in the region 
where the cut-off is O(Aqcd), is yet another approximation scheme which elucidates some 
of the observed features of hadronic physics. Strong coupling methods are usually employed 
in the context of gauge theories formulated on a lattice. In the strong coupling region, one 
can show that quarks are confined, chiral symmetry is spontaneously broken, the lightest 
hadrons have the correct quantum numbers, and so on. These features are expected to 
survive in the weak coupling region (i.e. as one proceeds toward the continuum limit by 
increasing the cut-off), since there is no phase transition separating the two regions. 
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Let us illustrate the idea behind the strong coupling expansion in the context of pure 
gauge theory on a Euclidean 4-dim hypercubic lattice with spacing a (see for instance [21]). 
The fundamental variables of this theory are the parallel transporters U {x,fi), associated 
with the oriented link from site x to site x + afi (see Appendix A for details of the notation). 
The Wilson action for this theory is defined in terms of the oriented plaquette variables 



-11 



x;fiu 



where the plaquette variables are defined in terms of the parallel transporters as 

Ux;^j.u = U\x, ijl)U\x + an, v)!! {x + au, fi)U {x, u) 
and g = g{a) is the lattice gauge coupling constant. The action Sw is often abbreviated as 

Nc 



S' 



w 



Re tr [Up - 1] 



where the sum includes every plaquette with both orientations. This is the simplest gauge 
invariant lattice action which reduces to the Yang-Mills action in the continuum limit 
a 0. 

To demonstrate confinement, one needs to show that it requires infinite energy to 
separate a quark-antiquark pair by an infinite distance. The expectation value of the 
Wilson loop, W (C) = (trf/(C)), probes the interaction potential between a static quark- 
antiquark pair. Consider C to be a large fiat rectangular loop with sides Ra and Ta. In a 
linearly confining theory, W (C) obeys the area law 

W (C) ~ exp {-aRTa'^) , 

with a non-zero string tension a. If the theory does not confine, then the Wilson loop 
would follow a perimeter law 

W (C) ~ exp [-/X {R + T) a] . 

At the leading order in the strong coupling expansion, one can easily show that Wilson 
loops obey area law as follows. The expectation value of the Wilson loop is 



WiC) 



1 



Y\ dU (x, /i) 



tr 



pX{u{i) 



exp 



2g^ 



where P denotes the orientation of loop C and 



Z' 



X{dU{x,^^) 



exp 



2g^ 
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Expanding the exponential in a Taylor series, we get 



exp 



Nc 



2^ 



mi ,m,2,... 



The group integration over [dll {x, fi)] produces a non-zero result only when all the factors 
of U {x, ji) for that link combine into a gauge singlet. The simplest instance is when one 
factor of U is combined with another factor of U\ 



JdU^l, jdU UijUl = ^J,k5u 



The leading contribution to W (C) thus arises when the smallest number of plaquettes 
(with matching orientations) tile the loop C. On explicitly carrying out the integration 
over the link variables, we get 

This shows that at strong coupling Wilson loops obey area law with string tension 

In (2^2) 



a2 



for Nr>3 



For Nc — 2, the orientation of the plaquette does not matter in group integration, and so 
each plaquette is associated with a factor oi This yields the string tension 

Subleading corrections to this strong coupling limit are conveniently obtained by expanding 
the exponential of the gauge action in terms of the characters of the gauge group. Such 
character expansions effectively sum a subset of graphs, and considerably simplify the 
computation. 

In contrast to the weak coupling expansions which are asymptotic in nature, strong cou- 
pling expansions have a non-zero radius of convergence. Within its domain of convergence, 
one can rigorously show: 

• Static quark confinement: Wilson loops in the fundamental representation obey area 
law at strong coupling. 

• Existence of a mass gap: The plaquette-plaquette correlation function decays expo- 
nentially at strong coupling. 

• Possible phase transitions: Power series expansions (in inverse powers of the lattice 
gauge coupling) for quantities like internal energy, mass gap, string tension etc. can 
be used to locate phase transitions. 



CHAPTER 2. QCD AT LARGE Nc AND STRONG TRANSVERSE GA UGE CO UPLING15 

When dynamical quarks are added to the gauge theory, calculations get technically more 
involved. Firstly, it becomes difficult to carry out strong coupling expansions to high 
orders. Secondly, phenomena associated with chiral symmetry breaking can only be seen 
if the quark part of the theory is taken into account exactly. This is due to the fact that 
chiral symmetry breaking produces pions as pseudo-Goldstone bosons, which are very light 
for small bare quark masses and hence give rise to long distance correlations. 

The results for the meson spectrum, at the leading order in strong coupling expansion, 
are described in Appendix C. They correspond to propagation of quark-antiquark pairs 
without any separation at all, and show that there is only one meson state for every spin- 
parity quantum number (i.e. all radial excitations are pushed away to infinite energy). 

2.3 Combined large- A^c and strong transverse 
gauge coupling limits 

We argued in Section 2.1 that large- A^c QCD in a good phenomenological approximation 
to many aspects of hadronic physics. This was deduced by power counting, combining 
the assumption of colour confinement for large- iVc QCD with the fact that at large Nc 
only planar diagrams with quarks on boundaries survive. In particular, to leading order 
in l/Nc, meson scattering amplitudes are given by sums of tree diagrams with exchange of 
physical mesons. The sum is always an infinite sum, because there are always an infinite 
number of mesons that could be exchanged in any given channel. This feature has a strong 
resemblance to the successful "Regge phenomenology" description of strong interactions. 
In Regge phenomenology, strong interactions are interpreted as an infinite sum of hadron- 
exchange tree diagrams, with the infinite number of hadrons lying on approximately linear 
trajectories in the — J plane. Thus it seems likely that the \arge-Nc limit of QCD 
should be an important cornerstone in the eventual derivation of Regge phenomenology 
[8n| . It is also believed that the straightness of Regge trajectories is related, by analyticity, 
to the narrowness of hadron resonances. If there is any approximation to QCD in which 
the trajectories are linear, this must be an approximation in which the resonances are 
narrow. Again it is an attractive feature of large- A'c expansion that the hadron resonances 
are narrow, with widths O (l/Nc). This offers hope that at A'c = cxd the Regge trajectories 
are linear. 

With the above in mind, the problem that needs to be addressed is how to sum the 
planar diagrams of large- Ac QCD. In order to tackle a quantum field theory, with its 
continuously infinite number of degrees of freedom, one has to impose kinematical cut- 
offs (i.e. a suitable regularization procedure). The dynamics of the degrees of freedom 
eliminated by the cut-off is accounted for by 

• renormalization of the existing couplings in the theory, and 

• introduction of new counter-terms with cut-off dependent couplings, that are consis- 
tent with the symmetries of the theory left unbroken by the selected cut-off. 
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The renormalized couplings and the induced counter-terms remove the cut-off dependence 
from physical observables and help restore continuum physics. The number and symmetries 
of the counter-terms depend on the symmetry properties of the cut-off. For example, in 
the weak coupling perturbative treatment of gauge theories, if one adopts a cut-off which 
maintains manifest residual gauge invariance (some gauge choice is essential to make the 
perturbative scheme well defined; usually a Lorentz covariant gauge is chosen) and Lorentz 
invariance at all intermediate stages of the calculation, then the counter-terms also respect 
the very same symmetries. 

In a confining theory such as QCD, gauge fields have large variations, and a non- 
perturbative cut-off should be such that it is manifestly gauge invariant. This criterion in 
the choice of cut-off is far more important than the criterion of manifest Lorentz invariance, 
because such a starting point is closer to the real world of strong interaction physics. (Of 
course, quantitative description of strong interaction physics will be attained only after 
taking in to account all the cut-off dependent counter-terms.) Such a cut-off is provided 
by formulating the gauge theory on a lattice. The discussion in Section 2.2 shows that for 
very coarse lattice spacing, the gauge theory produces linear confinement, i.e. large Wilson 
loops obey the area law. Numerical simulations for QCD, at different values of the lattice 
spacing, show that this behaviour persists all the way to very fine lattices with no phase 
transitions in the intervening region. 

An extremely coarse 4-dim lattice cut-off, however, has the disadvantage of drastically 
breaking Poincare invariance (actually 4-dim Euclidean invariance) all the way down to 4- 
dim hypercubic group. Though confinement is manifest, glueballs become infinitely massive 
and are pushed out of the spectrum. In the hadronic sector, only the lowest lying states 
survive in each quantum number sector. This is explicitly demonstrated in Appendix C, 
which analyses SU {Nc — > oo) lattice gauge theory at strong coupling with naive fermions in 
the fundamental representation. It has not been possible to perform analytical calculations 
on less coarse lattices; only numerical results have been obtained using powerful computers. 
Can we do something better without spoiling the analytical simplicity of working with a 
coarse lattice, and yet retain the higher energy excitations? Appendix B demonstrates that 
the 't Hooft model shares two key phenomenological features of strong interaction physics, 
namely, confinement and chiral symmetry breaking. Is there any way of incorporating 
these features in a systematic framework when studying (3 + l)-dim QCD? The answer to 
these two questions leads us to formulating QCD on a transverse lattice. 

To explicitly incorporate the 't Hooft model in the solution to large- Ai'c (3+l)-dim QCD, 
one must keep the temporal and one spatial dimension continuous. Crucial to the explicit 
solution of the 't Hooft model is the choice of the light-front = (x" + x^) / \/2 = as 
the quantization surface, and P" = (P^ — P^) / \/2 as the 'time' (a;+) evolution operator. 
Since P° > |P^|, it follows that P^ > 0. On the light-front, the 'energy' dispersion relation 
for a free particle of mass m is given by 



where p"*" is the longitudinal momentum and p± is the transverse momentum. Two impor- 
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tant observations follow from this dispersion relation: 

• = states have infinite energy unless p± = and m = 0. These 'infrared' diver- 
gences associated with longitudinal zero modes are special to field theories analysed 
on the light-front. These are separate from and in addition to the infrared problems 
encountered in a covariant formulation. The infrared divergences of the latter type 
show up in the light-front formulation as divergences associated with small p± states. 

• There are ultraviolet divergences associated with large p± states. 

One needs to regulate these divergences before any concrete calculation is attempted. The 
large p± divergences can be regulated by imposing a lattice cut-off a± in the transverse 
dimensions. We will regulate the longitudinal infrared divergences by explicitly eliminating 
the longitudinal zero mode using the principal value prescription. 

There are distinct advantages associated with this choice of cut-offs. 

• The total longitudinal momentum P+ of a state is just the sum of the longitudinal 
momenta of its constituents, P+ = J2i^t- It follows that the vacuum of the cut-off 
theory (for which P"*" = 0) contains no constituents, i.e. the vacuum is trivial. 

• In the = gauge, A~ is non-dynamical and plays the role of providing a linear 
confining interaction between quark-antiquark pairs in the longitudinal direc- 
tion. The infrared singularity of this confining potential is regularized by the principal 
value prescription. 

• After fixing the A'^ = gauge, one is still left with residual gauge degrees of freedom 
corresponding to a;~ -independent gauge transformations. These can be manifestly 
maintained on a transverse lattice. 

• Working with a very coarse transverse lattice makes confinement in the transverse 
directions automatic. Wilson loops in the transverse plane then obey area law. 

In our study of transverse lattice QCD, we restrict ourselves to the \aige-Nc and strong 
transverse gauge coupling limits [22]. This choice is dictated by the fact that only in these 
limits we can analytically solve the theory in a closed form. Consequently, our approach 
has following limitations: 

• Although linear confinement is built in to the theory, the dynamical mechanism 
behind it is different in longitudinal and transverse directions. Confinement arises in 
the longitudinal direction because the Coulomb potential in (1 + l)-dim is linear for 
any value of the gauge coupling, while it arises in the transverse direction because 
of the strong gauge coupling (i.e. completely disordered gluon fields). Thus perfect 
rotational symmetry is not expected in our results. 
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• Ghions have no transverse propagating degrees of freedom, so they can be completely 
integrated out of the theory. The glueballs are infinitely massive, just as in 4-dim 
strong coupling lattice QCD. 

• The spectrum consists exclusively of towers of meson states, and the large- A'^c limit 
brings out the constituent picture of mesons in an explicit way. But the large-A^,; 
limit also makes the baryons infinitely heavy. Although the baryons can be brought 
back in to the theory as semi-classical soliton states, that requires another set of field 
theoretical techniques. The case of baryons is not dealt with in this thesis. 

• Renormalization group scaling does not hold in the strong coupling lattice theory, 
and hence our results for dimensionless quantities are not expected to be the same 
as in the contimmm theory. The scaling behaviour can be improved by calculating 
higher order terms in the strong coupling expansion or by adding counter-terms to 
the simplest lattice action. We have not done that, and so our results can only be 
viewed as a model of real QCD. We hope that by having two dimensions already in 
the continuum, our results would not be too far off the real world, and our model 
would have practical use. 

• With the principal value prescription, the infrared cut-off introduced to eliminate 
the longitudinal momentum zero mode disappears from the final results. On the 
other hand, the transverse lattice ultraviolet cut-off that eliminates the transverse 
momentum modes with \p±\ > T:/a± explicitly appears in the final results. So our 
results can be interpreted only in the sense of an effective theory, i.e. valid at a 
scale below the cut-off scale where the contribution of the neglected counter-terms is 
suppressed. 

• Dynamical fermions on the lattice have the well-known confiict between chiral sym- 
metry and species doubling. For computational convenience, we have studied fermion 
actions with only the nearest-neighbour hopping term. The two popular choices in 
this case are: 

— Naive fermions: These retain exact chiral symmetry, but lead to multiple species 
of fermions — one for every corner of the Brillouin zone. Exact spin-diagonalization 
converts naive fermions to staggered fermions and reduces the number of fermion 
species, but that cannot eliminate the doublers completely. In the large-iVg limit, 
however, fermion doubling is not a serious problem because the fermion deter- 
minant does not contribute. 

— Wilson fermions: These avoid fermion doubling, but explicitly break chiral sym- 
metry. The "chiral limit" is defined as the location where pions become mass- 
less, but even at this location pions do not have all the properties of Goldstone 
bosons. 

Despite these limitations, our approach is tailor-made for the study of deep inelastic 
hadronic physics. Experimental results for deep inelastic structure functions are expressed 



CHAPTER 2. QCD AT LARGE Nc AND STRONG TRANSVERSE GA UGE CO UPLING19 

in terms of light-front variables, and their scaling behaviour implies that transverse dy- 
namics contributes only at the subleading order. So even if the transverse directions are 
severely distorted using a coarse lattice, the leading scaling behaviour can remain largely 
undisturbed. Note that for any 2-point or 3-point Green's function (e.g. propagators and 
form factors), it is always possible to choose a reference frame where the transverse mo- 
menta of the external states vanish. We make this choice in our analysis to express our 
results in the continuum language in a straightforward manner. 

2.4 SU {Nc) transverse lattice gauge theory 

Elements from both continuum gauge theory (Yang and Mills) and lattice gauge theory 
(Wilson) are needed in constructing transverse lattice gauge theory. Such a theory was 
first formulated by Bardeen and Pearson [l]. 

2.4.1 Basic gauge invariant quantities 

We start by introducing pure SU (Nc) lattice gauge theory on a transverse lattice in four 
dimensions. The underlying space-time is the direct product of two-dimensional Minkowski 
space and planar square lattice (with lattice spacing a±). 

Let Cy^yj_-x,xj_ be a particular directed curve from {x,x±) to {y,y±), and let the parallel 
transporter for the gauge field along the curve C be U {Cy^y^-x,x^) £ SU {Nc). Under a 
gauge transformation, 

U {Cy^y^■x,X^_) U' iCy^y^.x,xJ = V (?/, y ±) U ( , j^^ ;^ ) l^"^ (X, X ±) , (2.l) 

where V {x,x±) ,V {y,y±) are arbitrary SU {N^) matrices. For an infinitesimal straight 
segment from {x,x±) to {x + dx,x±), 

U {Cx+dx,x^;x,xi^ = 1 ~ (3^5 2:_l) dx^ + ■ ■ ■ , (2-2) 

where (x, x^) is the continuum gauge field, i.e. an element of the Lie algebra of SU {Nc). 
The gauge transformation rule for A^(x, x^) follows from the gauge transformation rule 
for parallel transporters, eq. ()2.Hl . 

A{x.,Xi_) A' {x,Xi_) = V {x,x±) A^{x,x±)V~^ {x,x±) 

+ V{x,x^)-d^V-^{x,x^) (2.3) 

Given the gauge field {x,x±), the parallel transporters can be constructed as, 

U{C) = Vexpl^-J^ A^{c{t),x^)—dt 
= Vexp!^-j A^{x,x^)dx'''^ , 
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where c^{t) parametrizes the curve C with c'^(O) = x'^ and c^(l) = y'^, and the symbol V 
denotes path ordering with respect to the variable t. 

The simplest gauge invariant quantity one can construct using the parallel transporters 
is the trace of the parallel transporter around a closed curve, also known as the Wilson 
loop variable. For a closed curve Co based at point {x,x±), 



ti U (Co) = tr 



(2.5) 



'Co 

When Co is an infinitesimal rectangle based at {x,x±), with sides dx'^ and dy^, 



U (Co) U {Cx,x^]X+dy.x^^ U {Cx+dy,x^\x+dx+dy,x^^ 

U (^Cx-\-dx+dy,x±;x+dx,x±^ U (^Cx~{-dx,x±;xx±) (2.6) 

= 1 - {df,A^ {x, x±) - d^Af, (x, xa_) + [A^ (x, xa_),A^ (x, x_l)]_ } dx^'dy" H 

The quantity in the curly brackets in the last line of eq. p.6|l is the field strength of the 
continuum gauge field 

Ff^u {x, x±) = d^A^ (x, x_l) - di,A^ (x, x_l) + [A^ (x, x_l) , A^ (x, x_l)]_ • (2.7) 



It follows from eq. ()2.1|l and eq. (|2.6|] that under a gauge transformation the field strength 
transforms as 

F^,u (x, x_l) F'^^ (x, x_l) = V (x, x_l) Ff,y (x, x_l) V~'^ (x, x_l) . (2.8) 

So the simplest local gauge invariant quantity in the continuum is 

iilF^"" {x^Xi^F"" {x,Xi_)]. (2.9) 

Let us now consider the finite straight segment from (x, x_l) to (x, x_l + a±n), and the 
corresponding parallel transporter U {Cx^x^+a^n;x,x^)- We abbreviate this transverse link 
variable as f/(x, xx,n). It gauge transforms as 

U (x, x_L, ra) — > U' (x, x_L, n) = V (x, x_l + a±^n) U (x, x_l, n) (x, x_l) . (2-10) 

The parallel transporter around the smallest lattice loop (i.e. the elementary plaquette), 

Umn {x, x_l) = f/^(x, x_L, m)?7^(x, Xj_ + a^_m, n)U{x, x± + a±n, m)U (x, x±,n) , (2.11) 

gauge transforms as 

Umn (X, X_l) f/^„ (X, X±) = V (x, X±) Umn (x, X±) V'^ (x, X_l) • (2.12) 

So the simplest gauge invariant quantity one can construct on a planar square lattice is 

tr \Umn 

{x,x^)]. (2.13) 
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In a transverse lattice space-time, we also need to consider "mixed" loops which traverse 
partly in the continuum and partly on the lattice. For a mixed loop C„ based at {x,x±), 
with sides dx and a±n, 

U (Cm) = {x, X±, n) U {Cx,x^+a^n\x+dx,x^+a^n) U {x -\- dx, X_L, Tl) U [C x+dx ,x ^]X ,x ^) 

= 1 + U'' (x; x±,n) {dpU {x; x±,n) + Ap (x, x± + a±n) U (x, x±,n) (2-14) 

—U {x, x±, n) Ap (x, x±)} dx^ + ■ ■ ■ 

The quantity in the curly brackets in the last line of eq. ()2.14p is the mixed gauge covariant 
derivative of the transverse link variable 

DpJJ (x, x±,n) = dfJJ (x, x_L, n) + Ap (x, x± + a_Ln) U (x, x_l, n) — U (x; x±, n) Ap (x, x_l) . 

(2.15) 

From the gauge transformation properties of parallel transporters, eq. ()2.3|) and eq. p.lOp . 
one finds that 

) — > V (x, x_L + a_\_n) D^U (x, x±,n)V ^ (x, x_l) , (2-16) 
and hence the simplest two-derivative gauge invariant quantity of "mixed" type is 

tr (^[D^U (x; x^, n)f D^U (x; x±, n)) . (2.17) 

Now let us introduce fermions, in the fundamental representation of SU (A^c), in to the 
theory. Under a gauge transformation V (x, x_l), the fermion fields \l/ and ^' transform as 

\E'(x, x_l) (x, x_l) = (x, x_l) \E' (x, x_l) 

^ (x, x_l) ^ (x, x_l) = {x,x±)V'^ {x,x±) . (2.18) 

The simplest gauge invariant fermion bilinear terms are therefore of the form 

,y±;x,x± )^(x,xx). (2.19) 

In a transverse lattice space-time, there are two distinct choices for the curve C in this 
term: 

1. C is an infinitesimal segment from {x — dx,x±) to (x, x_l). Expanding the fermion 
gauge invariant term, we get 

^ (x, x_l) U {Cx,x^;x-dx,xJ ^ (a; - dx, X_l) 
= '^{x,x±){l-Ap{x,x±)dx^'}{'^{x,x±)-^p^{x,x±)dx''} + ■■■ (2.20) 
= (x, x_l) ^ (x, x_l) - 1' (x, x_l) [dp + Ap (x, x_l)] (x, x_l) rfx^ H . 

The first term in last line of the above equation is the fermion mass term. In the 
second term, the combination Dp = dp + Ap is the gauge covariant derivative operator 
for fermions. Gauge covariant derivatives of the fermion field transform in a manner 
identical to the fermion field, 

Dp"^ (x, x±) ^ D'^^' (x, xx) = V (x, x±) Dp^ (x, x^) • (2.21) 
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2. C is the finite straight segment from {x,x±) to {x,x± +a±n). The corresponding 
gauge invariant fermion bilinear term is the transverse hopping term, 

1' {x, x± + a±n) U {x, x±, n) ^ (a;, • (2.22) 



In summary, the maximally localized gauge invariant terms involving the gauge fields 
Afj, (x, x±), U (x, x±, n) and the fermion fields \i/ (x, x_l), \I/ (x, x±) are: 

• ti{F^''{x,XA_)FP''{x,xi_)) 

• ti{\^D^U{x,xx,,n)]^ Di,U{x,x^,n)^ 

• tr(f/^„(x,x_L)) 

• \i/(x, X_L)\i/(x, X_l) 

• ^ (x, x_L + a_Ln) t/ (x, x_L, ra) \i/ (x, x_l) 



2.4.2 Transverse lattice QCD action 

Using the above terms, and demanding two-dimensional Poincare invariance as well as pla- 
nar square lattice symmetry, we can write down a general action for the SU (Nc) transverse 
lattice gauge theory: 

S = ^y"t/2x^tr(F^'^(x,x^)F^,(x,x^)) 

X_\_ 

+ f^Xci fd^xJ2^^{[D^U{x,x^,n)]^D^U{x,x^,n)^ 

^ xj_,n 

N r 

+ tttV^" Retr[[/^„(x,xx) - 1] (2.23) 

+ ai (i^x^l^(x,x_L) [^t'^-D;, - m] ^ (x,x_l) 

+ J (fx \^ (x, x_L + a±n) (r — ^7") U (x, x_l, n) \l/ (x, x±) 

x^,n 

(x, x_l) (r + 27") f/^ (x, x_L, n) (x, x_l + a±_n)\ . 

Here we have adopted the following conventions: 

• The normalizations for the gauge and the fermion fields are chosen such that each 
term in the action is proportional to Nc- With this choice, the large-iVc limit amounts 
to finding the stationary point of the action, while holding g fixed [TT] . 
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• Anisotropy parameters A^, A^^, k have been introduced to take care of explicit break- 
ing of rotational symmetry. 

• The Wilson parameter r has been introduced in the fermion transverse hopping 
term. As aj^ — 0, a non-zero value for r adds an irrelevant operator to the fermion 
Lagrangian, a_L(£'n^)(-D"\I/). That removes the fermion doublers without changing 
the continuum limit, but breaks the chiral symmetry explicitly. As mentioned earlier, 
we study two special cases: 

— r = 0: Transverse lattice QCD with naive fermions. 

— r = 1: Transverse lattice QCD with Wilson fermions. 

In order to reveal the connection of the above action with continuum 4-dim SU (Nc) gauge 
theory, we take the naive continuum limit a± ^ 0. The transverse link variable can be 
parametrized as 

^2 

U {x, x±, n) = exp [-a^A^ (x, = 1 - a±^„ (x, + ^-l) H (2-24) 



In the pure gauge field transverse lattice terms, using the forward lattice derivative 

f ^ _ Am (x, xx. + a±n) - Am {x, x^) 

and the Baker-Campbell-Hausdorff formula 



(2.25) 



exp(a;) exp(y) = exp ( x + y + y]_ H ) , (2.26) 



one finds that 



D^U (x, a;_L, n) = -a± {F^n {x, x±)} + 0(ai) ,^ 
Umn{x,x±) ^ex.p{-a'iFmn{x,x_L) + 0{a\)} ' 



with 



F^n {x, x±) = d^An (x, Xa_) - A{^^ (x, Xa_) + [A^ (x, X_l) , An (x, Xa_)]_ 
Fmn {x, xa_) = A4A„ (x, X_l) - A^A^ (x, X_l) + [A^ {x, X_l) , An (x, X_l)]. 



(2.28) 



Thus the leading terms, as a±_ 0, coincide with the continuum Yang-Mills action, 
d^x ^ tr ([D^U {x, x^, n)]^ D^'U (x, , x^, n)^ 

xj_,n 

J d'x J2 [«ltr {F^n {x, x^) F^^ (x, x^)) + 0(ai)] , (2.29) 
d'^x ^ Re tr {x, x±) - 1] 

x± ,m,n 

/a e 



x±,m,n 



^tr (F-" (x, x^) Fmn (x, x^)) + 0(ai) 



(2.30) 
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In case of the fermionic transverse lattice terms, the a± ^ limit yields, 

- / d^x a_L {x, x± + a±n) (r — ^7'*) U {x, x±,n)'^ {x, x±) 

x±,n 

+^ (x, x±) (r + 27") f/''" {x, x±,n)^ {x, x± + a±n)] 
= d^xJ2(^l I 7" A;^* {x, x±) - {Al'^ {x, x±)) 7'^* {x, x±) } 

x±,n '- 

+i {x, x±) 'j^An {x, x±) ^ {x, x±)] (2.31) 
+ Jd^xY, «i 2:^) * (x, xx) - ^ (A^^ (x, x^)) (A^* {x, x^)) + O (ai) . 

x±,n '- 

Thus the terms proportional to r shift the fermion mass by 2Kr/a±, and add an irrelevant 
operator to the fermion action, while the other terms agree with the continuum Dirac 
action. 

2.4.3 Determination of parameters 

The preceding analysis shows that in the continuum limit, all the anisotropy coefficients 
should tend to unity, 

a±^0: Xci, Xii,K^ 1. (2.32) 
The transverse lattice action then reduces to the conventional continuum action 

jS" = Sq -\- Sp 

^ ^ / / d^x^[F^^{x,xj_)F^,{x,xj_) + F''^{x,xj_)F^,,{x,xj_) 

+F'^''{x,x^)Fmn{x,x±)] (2.33) 
+ J d^x J d^x± ^ {x, x±) [iY ((9^ + A^ {x, x^)) 

+i7'' {dn + A„ (x, x±)) - m] * (x, x±) . 

In the quantum theory, the gauge coupling g gets converted to the QCD scale Aqcd by 
dimensional transmutation, while the quark mass m is not a physical observable. They 
have to be determined using experimentally observed hadronic properties as input. 

When the transverse lattice cut-off is finite, quantum renormalization effects make 
the anisotropy coefficients different from unity, i.e. they become functions of a±. For 
small values of a±, these coefficients can be computed using perturbation theory, and they 
remain finite. In general, they must be determined non-perturbatively, by demanding that 
physical observables satisfy rotational symmetry as closely as possible. Convenient choices 
for enforcing rotational symmetry are isotropy of the energy-momentum dispersion relation, 
and degeneracy of multiple helicity states of hadrons with non-zero angular momentum. 
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Aqcd and quark masses have to be determined from hadronic properties as usual, and so 
more physical conditions have to be imposed to fix all the parameters of the theory. For 
Wilson fermions in particular, m contains a power-law divergence as a± — ^ 0, and therefore 
it always needs to be fixed using a non-perturbative criterion. The conventional choice is 
to define the "chiral limit" for Wilson fermions as the value of m (not necessarily zero) 
where the pseudoscalar meson mass vanishes. 

The strong transverse gauge coupling limit can be viewed as the extreme case of 
anisotropic renormalization group evolution, where the transverse lattice cut-off scale is 
lowered as far as possible. In this limit a± is O(Aqcd), and we expect the anisotropy 
coefficients for the gauge field terms (Ac/, Xu) to vanish. We also do not expect any phase 
transition to take place during the renormalization group evolution from small a± to its 
largest value. These expectations are based on the non-perturbative behaviour of 4-dim 
lattice QCD, and are supported by leading order perturbative calculations |5j that show 
logarithmic behaviour for Ac/, Xu. Thus the action we have studied in this thesis is 

+ a\ I (fx^"^ {x.XA^liYD^-m]'^ {x,Xi_) (2.34) 

+ — ^ / (i^x [\E' (x, xj_ + a_Ln) (r — 27") f/ (x, x_L, n) \E' (x, x_l) 

(x, x_l) (r + Z7") U'^ (x, x_L, n) \E' (x, x±_ + a±n)\ . 

Compared to the continuum action, eq. ()2.33|) . this action has only one more parameter, 
i.e. K. Since the powers of a± that accompany factors of and k are obvious, in order to 
simplify expressions, we choose the units such that a_L = 1. 



2.5 The generating functional 

The hadronic physics is buried in various quark-antiquark correlations functions. As ex- 
plained in previous sections, large-iVc and strong transverse gauge coupling expansions 
combine efficiently to produce reasonable mesonic structure in the leading order. This 
section is devoted to computing the generating functional of quark-antiquark correlation 
functions in these limits. 

Since we are primarily interested in the quark-antiquark correlation functions, we add 
a suitable source term to the action, eq. (l2.34|] . 

i = S + j (fx£y Jtf {x,xr,y.y±)^t {y.y±) ■ (2-35) 

The generating functional, Z [J], is given by 

Z[J]= I [DA] \DU] [Dm ■ DW^ exp [il + (gauge terms)] , (2.36) 
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where "gauge terms" stand for the gauge fixing terms as well as the corresponding ghost 
determinant. The time-ordered quark-antiquark correlation functions can be obtained by 
taking functional derivatives of W [J] = —ilnZ [J] with respect to J. For example, 



^ dJij {x,x±;y,y±) 



(2.37) 



j=o 



T (< (x, x^) *f (y, y^) ^ {z, z^) {w, w^] 

T (^," (x, x^) v^f {y, yS) ) <T {z, z^) *f {w, w^)) ) (2.38) 



d'^W [J] 

—i 



djff (x, x±; y, y±) dJ^i (z, z±; w, w±) 



So if one is interested in computing the quark propagator and the quark-antiquark scat- 
tering amplitude, then W [J] should be known to O (J^). Before presenting the technical 
details of the calculation, we outline the essential features of the computation of W [J]. 

1. In the strong transverse gauge coupling limit, the dynamics of different x^-hyperplanes 
does not completely decouple. Colour singlet hadron states can jump from one hy- 
perplane to the next, and the fermion transverse hopping term (with coupling k) 
continues to contribute in this limit. 

2. We work in the light front gauge = 0. There are no ghost terms in this gauge. 
Also the functional integral over A^ becomes purely Gaussian, and can be evaluated 
exactly. This integration shows that the sole outcome of the light-front longitudinal 
gauge field is to mediate a linear Coulomb interaction between quarks and antiquarks, 
when they are propagating in a given hyperplane. 

3. The integration over transverse gauge link variables can be performed exactly in the 
large- A^c limit. This results in non-local interactions between local fermion bilinears. 

4. The fermionic functional integration is carried out in the large-iV^ limit by rewriting 
the integral in terms of non-local bosonic variables. The lavge-Nc limit of the integral 
is then equivalent to the stationary point limit of the bosonized functional integral. 

Now we demonstrate this calculation in detail. By scaling all fields by their canonical 
dimensions using the lattice spacing a±, one can work exclusively with dimensionless quan- 
tities: 

^/i — ^m/'^-L ~ ^-L/ct± (2.39) 

^' (x',x'_l) = (a_L)^ ^ {x,x±) {x',x'j_) = {aj_)^'^ {x,x^) 

U' {x'; x'j_, n) = U [x] x^, n) A!^ {x' , x'^) = a±A^ {x, x±) (2.40) 
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Henceforth we assume that this has been done, and 
AT r 

XJ_ 



+ I d?x^'^ (x, x^) [ij^Di^ - m] * {x, x±) 

+ — / (fx {x, x± + n) {r — ^7") U {x, x±, n) ^ (x, Xi_) 

x±,n 

+^ {x, x±) (r + i7") {x, x±, n) * (x, x±_ + n)] 
The functional integral in the light-front gauge, A'^ = 0, is 

Z [ J] = j [DA-] [DU] [D-^ ■ D^ 



(2.41) 



X exp 



(2.42) 



X exp 



X exp 



+i / d^xd^^ J] (x, xx) J^'^ (x, XX ; y, y±) *f (y, y±) 
j[d'xJ2 {^i ^± + - ^7")"^ {x; x^, n) *f {x, x^) 

x±,n 

(x, xi_) (r + i7")"^ 4 {x; x^, n) (x, x^ + n^) } 
~\J (^-^-^SuSjkd'iA-{x,x±)^ 

Xj_ 

- (fxY^ {^i 7+*j (a:^, 3;^)) Ar. {xxj_) 

The Gaussian functional integral over the A~ variables produces a non-local interaction 
term, which is quadratic in fermion bilinears (we suppress the overall normalization con- 
stant that does not affect expectation values). 

Z [ J] = J [DU] [D^ ■ D^ 

i (fxY^ {x, x±) {ij^df, - m)"/ *f {x, x±) 

+i / (fx(fy ^ *"(a;,a;x) J^^(a;,xx;y,yx)*^(y,?/x) 

■ d'xd^y Y: h^r' Sx.y.h {x - y) (2.43) 



X exp 



X exp 



2N. 
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{x, x^) {y, y^) {y, y^) (x, x^] 



X exp 



Y j (fx ^ (x, x± + n) {r - 27")°^ Uij (x; x±, n) (x, x_l) 
(x, x±) (r + i7")"^ U], (x; x±, n) (x, x± + n) 



where h{x — y) satisfies d'^h{x — y) = 5^'^\x — y). The next step is to carry out the functional 
integral over the U variables. In the large- Ai'c limit [2] [T?^ . 



Z[J] = J [D'^H-DW] 



X exp 



d^a; J2 < (a;, a;^) (zt'^S^ - m)^^ (x, x^) 

Xj_ 

+i j d^xd^y ^ ^'"(x,x±) J?^^(x,x±;?/,?/±)^^(?/,?/x) 



X exp 



2N, 



d'xd'y J2 {l^y'{lYs^,y,hix-y) 



(2.44) 



x±,y± 



X exp 



< (x, xx) (y, ^ (y, y±) (x, xx) 
-NcJ d^x J2 j Jl + ^5 (a;; x^, n) - 1 

xj_,n I. V 



1 1 



where B (x; x_l, n) is a matrix in Dirac space with matrix elements 

B"'^ (x; x_L , n) = -"^^ (x, x_l + n) (x, x± + n) {r + ^7^) ^'^ 
• ¥J (x, x_l) (x, x_l) (r - 27^^)^'^ 



(2.45) 



To perform the functional integral over the fermion variables, we use a simple trick (see 
Appendix B, eq. (IB.8jl ). We first introduce auxiliary non-local bosonic variables to rewrite 
the four-fermion interaction term in terms of fermion bilinears, and then integrate out the 
fermion variables. 

Z[J] = J [Da ■ D\] exp [iV (a, A; J)] , (2.46) 

where 

a''''ix,xr,y,y±) = ^"ix,x^)<il^iy,y,_) . (2.47) 
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Then 



V [cr, A; J] = — /i [a] + i/2 [c] — J (fx Tr^ln [—iX{x, x±;x, x±)] ^ 

/ 



Kj y^) ^tj (^^ y^ y±) 



(2.48) 



+af {x, xr, y, y±) {S^^y^ {zYd, - m)^f 5' {x - y) + Jff {x, xr, y, y±)) 



with 



MB] 



2Nr 



d'xd^Y.hT'h^y' 



Sx^y^h {x - y) (x, x^; y, y±) a]^ {y, y±; x, x±) 



7/3 



^ Nc f d^x J2 ti-D (/2 [B {x; x±, n)] ) , 



1 1 



1 + ^B-l-ln - + -W1 + — . 



2 2 



Nl 



(2.49) 
(2.50) 

(2.51) 



The matrix elements of B expressed in terms of a are, 

(a;;x_L,n) = -a'^^ {x,x± + n;x,x± + n) [r + i'y'^Y^ 
■ cr]J {x, x±; X, x±) (r - i7^)^'^ . 



(2.52) 



In the large- A^c limit, [cr, A; J] is of 0{Nc). Therefore, the Nc ^ 00 limit amounts to 
evaluating the functional integral at its stationary point. The stationary point {a, A) is 
determined by 



dV 



dX^f {x,x±;y,y±) 
dV 



da"f {x,x±;y,y±) 



(J, A 



-a/3 



= A,7 y^ y^) + ^^±yi. (^t^^m - ^Ti - ?/) 



+Jtf {x,x^;y,yA_) - —6ij{-i+Y b'^Y - y)^x^y±(^ll{y,y±;x,x±) 



-iN,xl,j(9 [x - y) 5,, 



(2.54) 



• ^ I (r + i7^) o-fefe (x, x± - n; x± - n) (r - i-fl) [B (x; x± - n, n)] 

n 

+ - ^7n ) /2 [-S (2;; 2;±, n)] Ukk {x, x±_ n; x, x_l n) (r i7^) |^ = , 
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where 



The first stationary point equation ()2.53|) is trivially solved by 



-1 /3a 



X^^ {x,x±;y,y±) =i{a ^)'^." iy,y±; x, x±) . 
Substituting this in the second stationary point equation p.54|) . we obtain 



(2.55) 



(2.56) 



i ((^ ^Yii {y, y±; x±) + Sx^yj_ {il^d^ - m)"^ 5^ (x - y) 



af3 ^2 



+Jif {x,x±;y,y±) - j^6ij {-f^)"^ {l^Y"^ h{x - y)6^^y^all{y,yr,x,x±) 
-iNJij5'^ {x - y) 5^^y^ 

■ ^{{^ + nj) o'kk {x, xj_ - n; X, xj_ - n) {r - ij'^) [B (x; x± - n, n)] 



(2.57) 



+ {r- ill) f'2 [B x^.n)] a^k {x, x± + n; x, x± + n) {r + i'jl) 



13a 







Solution of this equation will yield a as a function of J. The large- Ai'c generating functional 
can then be succinctly written as 

Z [J] = exp [iV {a, A; J)] W[J]=V {a, A; j) = Kfr (^; J) , (2.58) 

where Kfr (o"; J) is obtained by substituting A in terms of a. Modulo J-independent terms, 

Ves (a; J) = -/i [a] + 1/2 [a] + i [ d\ ^ [Tr (In a) (x, xr, x, x^)] (2.59) 

£x£y alf (x, x±; y, y±) S^^y^ {i^d^ - m)^ 5^ {x - y) + J^f (x, x^; y, y±) . 

x±y± 

The generator of IPI vertex functions T [ip] (also known as the effective action) is the 
Legendre transform of Let the effective field ip be the variable conjugate to the 

external source J. Then 

a/3 . ^ SW [J] 

y^ij {x,x±;y,y±) - 



5Jif ix,x±;y,y±) 



(2.60) 



T[ip] = W[J]- (fx(fy Y] v^"^ (a:,x_L;?/,?/_L) (x,x_l;|/,|/±) 



(2.611 



+ I £xd^y Y "Pif i^^ y±) [^x^yj_ {il^d^ - m)°^ 5^ (x - y) 
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In the large- A^c limit, (p — a + O {l/Nc). It is also easy to show that 



where we have used the abbreviations, 
S^W [J] 



ij;kl 



{xx±,yy±;zz±,ww±) , 



(5J°^(a:, xr, y, y±)5J^f{z, zr, w, w±) 



= [r^'^] ^^^) 



Siplf{z, z±; w, w±)S(pfJln{u, u±; v, v±) 



(2.62) 
(2.63) 

(2.64) 
(2.65) 



Chapter 3 

The Chiral Condensate 



ofree 



3.1 Chiral symmetry 

Consider the transverse lattice action for free naive fermions 

/ Sx ^ [-i^ (x, x_L + a±_n) 7"\[^ (x, 

(x, x_l) 7""^!/ (x, x_L + a_Ln)] . 
On-shell fermion modes of this action satisfy the dispersion relation 



+ 



(3.1) 



E{p\p^)=± 



2 ^ 

{p^) + — sin^ (p_L ■ na_L) + 



(3.2) 



where E and (p^,p±) refer to the energy and three-momentum respectively. Since the 
transverse momenta n" range over the Brillouin zone, — — < n" < — , it follows that 



E{p\p^)=E{p\p^ + ql) 
where q]_ is one of the four transverse momentum vectors 

TT \ 



9l= (0,0), — ,0 , 0,— , — ,— 

I V«-L /V «-L/ V«-L a±/ 



TT \ 



(3.3) 



(3.4) 



The consequence of ()3.3|) is that in the naive transverse lattice fermions represent four 
fermion states per field component. This is the generic "fermion doubling" phenomenon 
on the lattice. In order to get a single fermion state per field component, one can add the 
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irrelevant Wilson term to the naive fermion action (cf. (l2.23jl ). 



+ 



/ (fx^^\^ (x, a;_L + ai_n) (r - Z7") ^ (x, x^) (3.5) 



+^ (x, x_l) (r + ^7") \1/ (x, x_L + a_Ln) 
With the r-dependent terms, the dispersion relation for on-shell fermion modes becomes, 



{p^Y + ~ si^^ [P-L ■ ""-^-l] + ( ^ cos [p_L ■ na_L] ) 

n -'- \ n / 



. (3.6) 



This result shows that the p± = mode at the centre of the Brillouin zone is massless at 
m = 2Kra'J^. Moreover, compared to this mode, the other p± = q']_ modes get an additive 
contribution n^2Kraj[^ to their masses, where is the number of components of equal 
to 7r/a±. The modes with 7^ thus become infinitely heavy in the a_L — > limit, leaving 
behind only one physical mode with finite energy at p± = 0. The usefulness of Wilson 
fermions comes to fore when there are interactions between fermions. In an interacting 
theory, the fermion doublers, if not eliminated, can affect the physical content of the theory 
in a non-trivial way. 

Returning to the case of naive lattice fermions, one can show that their group of chiral 
symmetries is considerably larger than that of the continuum theory. Let us find the 
symmetry groups of individual terms in the naive fermion transverse lattice QCD action. 
In terms of variables scaled using the lattice spacing a±, 



S 



NF 



+ 



/ (fx'^'^ (x, x_l) [i'l'^D^ - m] ^ (x, x_l) 
— J (fx [— (x, x_L + n) 7"[/ (x, x_L, n) \l/ (x, x±) 



(3.7) 



(x, x_l) 7"'[/^ (x, x_L, n) ^ (x, x± + n 
With x_L = (x^,x^), let us define new spinor variables x(a;, x_l) and xi^i^A 

^ {x,Xi_) = {-i'^f\-i^f\{x,Xi_) , 

'^{x,XA_)=x{x,XA_){-i^f {-i'^f . 

In terms of these new variables, the action takes the form. 



(3.8) 



S 



NF 



X, X_|_ , 



m 



X(x,x_lJ 



+ j d^x^^n [-ix (x, x±+n)U (x, x±,n)x {x, x± 

+ix (x, x_l) (x, x_L, n) X (x, x± + n 



(3.9) 
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where 



A„ (x^) ^ A„ (x^ x') = <! [ ^ ^ - 2 ' (3.10) 

(-1) : n = 3 . 

In this form, the naive fermion transverse lattice QCD action has been spin-diagonalized 
in the transverse spin space. Let a lattice site be defined as odd (even), if x"^ + x^ is odd 
(even). The fermion transverse hopping term has an exact global symmetry Uo{4:) (g) ^4(4), 
defined by the transformations 



X {x, x±) x' {x^, x±) = X {x, x±) [/J(^) 



for x± odd (even). (3-11) 



The fermion mass term is invariant only under the diagonal subgroup Ua{4:) of [/o(4)(g)[/c(4), 
given by Uo = U^. To see the symmetry group of the continuum fermion kinetic term, we 
choose a particular gamma matrix representation: 

' ^ ~ V la^ ' ~\ia^ ) ' ^ ~ \ -<y^ 

We also decompose the fermion spinors in two component notation: 

X=(^j)> X=(^P ^q), (3.13) 
In terms of 2-dim spinors and ^q, the naive fermion transverse lattice QCD action is 
Snf = j dFx^'^x{x,XA_) i{a^DQ + i(T'^Di) -{-If^^'''' rr^'^),{x,XA_) 

+ y d^x An {x±) [-i^x {x, x± + n)U {x, x±, n) {x, x±) (3.14) 

X,x±,n 

+i'^x {x, x±) {x, x±, n) *A {x, x± + n)] , 

where the index A takes two values, P and Q. It is now easy to see that the continuum 
fermion kinetic term is invariant under the [7^(2) (8) Ur{2) transformations: 

^x [eyi]){ie''yT^®l+ie\-T^®a\-lf^+''')]^^,^!x'. 

*A ^ ^A' [exp(-i^^- r'^ 1 + i^^-r^®a3 (-1)^2+^3^]^,^, (3.15) 

where r'^ {k — 0, 1, 2, , 3) generate the Lie algebra of the group U{2) acting on index A. 

In summary, the symmetries of individual terms in the naive fermion transverse lattice 
QCD action are: 

• Fermion transverse hopping term (with coefficient k): Uq (4) <S) Ue (4). 

• Fermion mass term (with coefficient m): ^7^,(4). 
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• Continuum fermion kinetic term: Ul (2) ® Ur (2). 

When m = 0, the combined chiral symmetry of the action is Ul (2) ^Ur (2). When m 7^ 0, 
only the diagonal part of Ul {2)®Ur (2), i.e. Uy (2), remains an exact symmetry. The non- 
diagonal part of the symmetry is broken explicitly by the mass term, and spontaneously 
by the strong interactions (as we will show in the following sections). 

Eq. ()3.15|) also demonstrates that the doubling problem for naive fermions can be made 
less severe by dropping the sum over A. This gives the staggered fermion transverse lattice 
QCD action, which has two component spinors, and which gives rise to two degenerate 
fermion flavours in the continuum limit. 



3.2 Quark propagator 

The quark propagator is obtained from the generating functional of connected quark- 
antiquark correlation functions using 



T(xl/r(x,xx)vl/f {y,y^ 



dW [J] 



dJff (a;,x_L; y,yA_) 



(3.16) 



j=o 



where W [J] = V^ff (a; J) is given by eq. ()2.59p . In terms of the solution to the stationary 
point equation ()2.57p . 



T (^^" (x, Xi_) ^'J (y, yx) j ^ = (fXo)°/ (x, xx; y, y±) , 
with ao = a(J = 0). Setting J = in the stationary point equation (I2.57|l . we get 

^ i^o^Yii ^' + ^^±y± (^7^^M - "^)? ^^^^ - y) 



(3.17) 



(7+)°'^ {-i^y^ h{x - y)5.^^y^ (cro)Ifc (z/, Z/±; a;, x^) 



x; x± — n,n 



+ {r- ill) f2 [Bo {x; xx,n)] {ao)kk i^, x± + n; x, x± + n) {r + 27^) 



I3a 



(3.18) 







where Bq = -B(cro)- We expect ctq to obey (1 + l)-dim Poincare invariance and lattice 
translational invariance, and so we parametrize it as 



(^o)r/ {x,xr,y,y±) 



1^ 



■ 6 



-ip-{y-x)+ipj_-{yj_-xj_) 



13a 



m — H {p,p±) + ie 



(3.19) 
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In the coincidence limit x ^ y, the quark propagator is then proportional to the identity 
matrix, 

(^oXf (x,xx;x,a;x) =iV,ao5"^. (3.20) 

This coincidence limit is singular, and must be regulated. The regulator must be gauge 
invariant, and preferably also mass independent. A reasonable choice is the split-point 
regularization used in studying the chiral condensate in the 't Hooft model 

Assuming that 5"o has been properly regulated, we substitute ()3.2()ll in ()2.52l) and (j2.55|] . 
and obtain 

Bo (x; x^,n) = N^a^ (l - r^) , (3.21) 



/2 [Bo] 



K 



2N^ (1 + + (1 - r2) K^a^ 



Gi (gp) 
N3 



(3.22) 



where Id is the identity matrix in Dirac space. Using these results, the stationary point 
equation (|3.18|) becomes 



j^S^j (7+)"'^ (7^)^^ h{x - y)6^^y^ icro)ll {y, yr, x, x^) 



(3.23) 



(x - y) 5.,,,r^4 (1 - r^) aod (ao) = . 

With the parametrization ()3.19|) . this stationary point equation can be converted to an 
integral equation for S (p, p_l). 



{2nf [{p 



7 



1 

— m — T, (q) + ie 



7 



(3.24) 



- 42 (1 - r^) (ao) 



This is the same equation as in the 't Hooft model (see (IB.22|l ). except that the quark mass 
has been additively renormalized by a constant that depends on the chiral condensate. It 
can be solved exactly as 



4(1- r^) zaoGi (^o) • 



Putting everything together, the quark propagator takes the form 



T[^, {x,x^)^^{y,y^: 



(2vr: 



x,x±;y,y±) 
1 



(3.25) 



(3.26) 



(m-4{l-r')zaoG^ (^o)) + + 



I3a 



In particular, for Wilson fermions (r = 1), the additive correction to the quark mass 
vanishes. The quark propagator then reduces to that for the 't Hooft model (see (|B.25jl ). 
This feature is a consequence of the fact that the Dirac space structure for the Wilson 
fermion hopping term has projection operator form. 
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3.3 Chiral condensate 

Spontaneous chiral symmetry breaking is an important phenomenological property of the 
strong interactions (see for instance Approximate chiral symmetry, SU ® SU r(2) , 

arises in QCD, because the two light quarks, u and d, have masses much smaller than Aqcd- 
(Inclusion of the somewhat heavier s quark, extends the approximate chiral symmetry to 
SU L^'i) ® SU r{?i) .) Unbroken chiral symmetry is undesirable, because it predicts nearly de- 
generate parity doublets for hadrons — a property that is in contradiction with experimental 
facts. The conflict is avoided by realizing that the QCD vacuum spontaneously breaks the 
S'[/l(2) ® SUr{2) symmetry to its diagonal subgroup S'?7v'(2), i.e. in the same direction 
as the explicit breaking produced by the quark mass term. This pattern of spontaneous 
symmetry breaking predicts the existence of nearly massless pseudo-Goldstone boson with 
the same quantum numbers as the broken axial symmetry generators. In fact, the lightest 
hadrons are pions with precisely these quantum numbers, and one identifies the pions as 
the pseudo-Goldstone bosons. A non-zero vacuum expectation value of as m — 0, is 
an indicator of spontaneous chiral symmetry breaking. Demonstration of this behaviour 
in QCD involves all the complications of strong interaction dynamics, and we now analyse 
this feature in strong transverse coupling QCD. 



3.3.1 Naive fermions 



We have described the chiral symmetry properties of transverse lattice QCD in section 
(13.111 . Only for fermions with r = 0, chiral symmetry can be spontaneously broken. We 
thus have to demonstrate that that the chiral condensate (^'^'^ is non-zero in the chiral 
limit for naive fermions. From eq. ()3.26p . it follows that (which is the chiral condensate 
up to a proportionality factor) satisfies the recursive relation. 



(2^ 



tlr 



where, 



i» - (m - 4zaoGi (ao)) + ^ + le 



K 



(3.27) 



(3.28) 



The integral over p can be easily evaluated, with the principal value prescription [TT] . 



zero 



<Pp 

m - 



m 



Aiaod ((To) 



2p+p- - (m - AtdoGi iao)y + ^ + 



—I 



Aiaod ((Jo) 



47r 



dp~^ 
p+ 



(3.29) 



The remaining integral over p^ is logarithmically divergent. This divergence is along the 
light-front, and needs careful regularization respecting Poincare and gauge symmetries. 
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In the limit m — > 0, and with the introduction of a regularization kernel K{p'^,(Tq), 
eq. ()3.29|) takes the form 



In case of the 't Hooft model, Burkardt et al. have evaluated the chiral condensate of this 
form, using split-point regularization and operator formulation with the result (IB.51jl . 
A similar procedure would give a non-zero chiral condensate in our case too, but we have 
not been able to evaluate it in a compact form. 

3.3.2 Wilson fermions 

For Wilson fermions, chiral symmetry is explicitly broken. The condensate (\I'\I') therefore 
does not have direct phenomenological relevance, but we can still calculate it. For r = 1, 
the quark propagator is identical in structure to that for the (1 + l)-dim large-A^'c QCD, 
except for the dimensionality of Dirac space. Therefore, when the same regularization is 
used in (1 + l)-dim and (3 + l)-dim theories, 




(3.30) 



(3.31) 



Chapter 4 

The Meson Spectrum 



4.1 Quark- ant iquark scattering Green's function 

The time-ordered four-point quark-antiquark Green's function can be obtained from the 
generating functionals as 



T (x, x^) xl/j (y, y^)) ) (T {% (z, z^) vl/f w^)) ) 
d^W [J] 



djf {x, x±; y, y±) dJl' {z, z±; w, w±) 



7<5 



,yy±', zz_i_,ww±) 



J=0 



Using the effective action ()2.6H) . we find 

r(2)| 



S + K 



ix^±,yy±; zz^,ww±) = -i ((To (y, y±; z, z^) {a^ ^)^" {w, w±; x, x±) , (4.3) 



(4.1) 



(4.2) 



d'{-hM + ih[v]) 



difif (x, x±; y, y^) d^li {z, z^, w, w^) 



■6ij6kih {x - y) 5x'xj/x'^^ i.^ - ^) ^xxt«x^^ {y - ^) ^y±z± {l^T (7^) 



-i^Sij6ki6^ {x - y) 5'^ {x - z) 5^ (x - w) 5^^yJ^^^^ 



(4.4) 
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Here Gi and G2 are functions of cxo given by 

,2 



Gi (ao) 
G2 (ao) 



AC 



2 1 



1 



1^1 — r^) k^cTq 



4v/l + (1 - r2) /«2a2 + ^1 + (1 - r2) fi:2^2y 
The full quark-antiquark scattering Green's function is 



(4.5) 
(4.6) 



T < (X, Xx) Vl/f (y, y^) % {z, Z^) [w.W^] 



(4.7) 



(cro)r/ (a;, xr, y, y±) {(Jo)11 {z, z±; w, w±) + i[S + K] 



\7'5 



ij-ki {xx±,yy±;zz±,ww±) 



The meson spectrum is determined by the poles of this Green's function, while the corre- 
sponding residues are related to the meson wavefunctions. Clearly, the singularities of the 
quark-antiquark scattering Green's function are contained in the zeroes of [5 + 7^]. The 
homogenous equation which determines these zeroes and the corresponding eigenfunctions 
is the Bethe-Salpeter equation. 



4.2 Bethe-Salpeter equation 

In a symbolic form, the Bethe-Salpeter equation is 

[S + K]x = 0. 

In a more explicit notation, 

(fzcfw ^ ('S')^'^J^ {xx±, yy±; zz±,ww±) 



(4.8) 



+ i^TijM i^^±^yy±'^ zz^,ww±) xTi {zzx,ww±) = o . (4.9) 



7<5 



Let us transform this equation to momentum space. It follows from eq. ()3.26jl that 

{^0^)1^ ix,x±;y,y±) = % 



p-ip-(y-x)<: 
2 ^ "a^xJ/± 



■ 5. 



(2vr) 

- (m - 4i (1 - r^) d^Gi) 



Define the momentum transform of Xi!^ {x-, x±] y, y±) by 



xtf ix,x±]y,y±) 



W (2vr) 



o 

xtf {ki, ki±; k2, k2±) e~ 



9^ 
27^p^ 



/3a 



(4.10) 



■iki-x+iki±-x±—ik2-y+ik2±-y± 



4 Aij 



(4.11) 
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With ()4.3|] . (I4.4|] . (I4.1()|) and (I4.1H1 . the momentum space Bethe-Salpeter equation for 
meson states becomes, 



A - (mi - 4 (1 - r^) iaoGi) + 



9 _L -I Q(5 

2Trkt 



]k2 - (m2 -4(1- r^) iaoGi) - 



7/3 



xt if^2, hi] h, ki±) 



+ 



(2^(M 



.^(4) (^^ + - /i - h) — % 



{kt-ity 



—I 



2 (Gi + (1 - r^) ao^Gs) ( [r^^^^^^^ + (7")"^ (7")^^] cos {{k,^ + k^^) ■ n) 

n 

+r [^°' (7")^'' - (7")"' 5^^] sin ((fci^ + k2A_) ■ n) 



-4 (l-r^j'a^G^rt^'^} 



Xii {h, h±'-, h, h±) — . 



(4.12) 



Rewriting the above equation in an obvious matrix notation, we get 

^i + (mi-4(l-r2)^aoGi) + -^7 



Xii iP2,P2±;pi,Pi±) 



X / ^^^6^'Hpi+P2-h-k, 



- (mi - 4 (1 - r2) ^aoGi) V f + 

1 

. . 

JJ 



{p\-kty 



(27r) 

2(Gi+(l-r2) a^Ga) 

5Z y^ii *^^2, ^2±; fcl, fcl±) + 7"X5 (^2, ^2±; ^l, fcl±) 7" cos ((fci_L + fc2± 



X] [^5 (^2, ^2±; /ci, fci±) 7" - 7"xjj (^2, /i;2±; A;i, A;i±)] sin ((A;i_l + /i;2±) ■ n) 



-4 (1 - r^) V0G2X5 (^2, fc2±; ^1, fci± 
- j62 + (m2-4(l-r2)iaoGi; 



(4.13) 



27rpJ 



■7" 



p2 - (m2 - 4 (1 - r2) ^aoGi) V ^ + «e 
Without loss of generality, we scale the longitudinal momenta as 

pX ^pt = pX = kf = y = H = y2, 



(4.14) 



and choose to work in a reference frame where the total incoming and outgoing transverse 
momenta vanish, 

Pi±+P2± = 0. (4.15) 
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In such a case, except for the external quark propagators and the eigenfunction x, all 
other factors in eq. ()4.13p are independent of the "-" and "±" momentum components. 
This instantaneous nature of the Bethe-Salpeter kernel makes it convenient to project the 
Bethe-Salpeter equation on the light-front by |2] 



$ (x) 



dPi dpl 



(27r)2 (27r)^ 



6 



Pi -P2 j ^^^^ iPl±+P2±) 



xl{P2A-x,P2±;Pi,x,pi±). (4.16) 
Performing a simple contour integral over Pi , we obtain the homogeneous equation 



_ (mi - 4 (1 - r^) ta,G,f - ^ _ jm, - 4 (1 - r') i~a,G,f - 



X 



1 — X 



^(x) 



2x(l — X 

dy \ 



(mi - 4 (1 - r^) zaoGif-^ 2..^r\^^"^^' 
+ X7 + (mi - 4 (1 - r ) la^Gx) + 



2x 



27[X 



2{Gi+{l-r')alG2) [2r2$(i/) + 7"$(z/)7"l - ^ {l - r')' a^G^m 
M2 (,„,_4(l-r2)^aoGi)'-2! 



2 2x 

-(1 - x)7" + (ma - 4 (1 - r^) iaoGi) 



27r(l -x) 

Solutions to this equation provide the meson spectrum and light-front wavefunctions. 



(4.1^ 



4.3 Meson properties 

The integral equation ()4.17p is quite complicated, and it has to be solved numerically in 
general. But by comparing it to the corresponding equation for the 't Hooft model, (jB.37j) . 
we observe that four changes have taken place in going from the (1 + l)-dim theory to our 
(3 + l)-dim case: 

• The meson wavefunction has become a matrix in Dirac space. Physical meson states 
labeled by spin-parity quantum numbers have to be obtained by diagonalization of 
the matrix integral equation. 

• The quark mass has been additively renormalized due to hopping in the transverse 
directions (tadpole diagrams), m — m = m — 4(1 — r'^)id-oGi. 
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• The transverse lattice dynamics has shifted the interaction kernel by a term, indepen- 
dent of the longitudinal momentum fraction x but dependent on the spin-parity of 
the meson. This wavefunction at the origin effect is a result of the strong transverse 
gauge coupling which does not allow the quark and the antiquark to separate. 

• The 5-function constraint in the transverse directions makes the meson orbital angu- 
lar momentum vanish. Then Jz = Sz, and with spin-half quarks, the meson he- 
licities are restricted to 0, ±1. Allowed spin-parity quantum numbers for the mesons 
are therefore = 0^, 1^. 

These changes are simple enough, so that even without explicitly solving (I4.17jl . we can 
infer some of the meson properties, based on the known results of the 't Hooft model 
and noting that extrapolation in k from zero to non-zero value is smooth. We find it 
convenient to describe the meson wavefunctions in the basis that is a direct product of the 
usual Clifford bases in the continuum and the lattice directions: 

$ = 5^ $C;Lr^'^ , (4.18) 
C,L 

Tm = Tc0Tl = |i,7+,7",^[7+,7"]} ® { 1,7^,7', ^[7',7']} • (4.19) 

In this basis, 7"'^c',i7" ^ ^c,l for each value of L. The sixteen components of ()4.17p . 
therefore, separate in to four blocks (labeled by the value of L) of four components each. 
In addition, parity remains an exact symmetry in each block. 
It is then easy to deduce the following features: 

• The singular part of the interaction kernel is the same as in the 't Hooft model. It is 
this part which determines the behaviour of finite norm solutions at the boundaries 
at X = and x = 1 [26l, and only the components contribute to it. We thus 
expect $-;L to vanish at the boundaries as and (1 — x)^^, where 

~ 2 ~ 2 

npi cot (7r/3i) = 1 ^ , vr/^a cot (vr/^a) = 1 ^ • (4.20) 

9 9 

• The light-front wavefunctions are gauge invariant by construction, and they are re- 
stricted to the finite box x G [0, 1]. The finite box size guarantees that the spectrum 
of is purely discrete In each spin-parity sector, meson states can thus be 
labeled by a radial excitation quantum number n = 1, 2, 3, . . .. 

• For large n, the behaviour of eigenvalues and eigenfunctions of the integral equation 
is governed only by the singular part of the interaction kernel. Then as in the 't Hooft 
model [2E], 

~ irg'^n 



<l>1.r(x) ~ V2 



sm nnx 



n>l. (4.21) 
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The non-singular spin-parity dependent part of the interaction kernel, arising from 
the transverse lattice dynamics, only provides an n-independent zero-point energy 
shift. We thus expect parallel infinite towers of meson states in the — n plane, 
labeled by spin-parity quantum numbers. 

• Under the exchange of the quark and the antiquark, the meson wavefunction trans- 
forms as Eqg^(x; 1711,1712) = C$^(l — x;m2,mi)C~^, where C is the charge conju- 
gation operator. Eq. ()4.17|l is invariant under this exchange operation (it is easier to 
see this before performing the contour integral over p^). Although parity is an exact 
symmetry of our formalism, it is not manifest, because the light-front is not invariant 
under parity transformation. The {x 1 — x) part of Egg can be associated with 
parity, however, which implies that meson states in each tower alternate in parity as 
in the 't Hooft model. Since is symmetric, and the quark and the antiquark 
have opposite intrinsic parities, the lowest meson state in each tower has negative 
parity and P = (—1)". 

• The parity symmetry is exact in every block labeled by L, so the eight possible spin- 
parity quantum numbers are distributed in to the four blocks as two states of opposite 
parity in each block. Explicitly, in the conventional notation, {7r,ai(0)} G ^C;nin2: 
{p{n),ai{n)} G $c;n and {p(0),a} G 

• Another consequence of the parity symmetry is that the equations for the components 
$C;ni and $c;n2 are degenerate. This feature implies identical masses for meson states 
with helicity= ±, which correspond to the transverse lattice combinations (ni ±^2) 
(e.g. in case of vector and axial mesons). There are thus three independent towers 
of meson states. Note that helicity can be defined only modulo-4 on a transverse 
square lattice, and explicit breaking of rotational symmetry makes the equations for 
the helicity= state and the helicity= ± states non-degenerate. 

• Fitting the meson spectrum to experimental results requires determination of the 
parameters and k, in addition to the quark masses. It would be convenient to fix 

using the slope of the infinite towers of states, and k by demanding (to the extent 
possible) degeneracy of helicity= 0, ± states. 

Now we consider the cases of naive and Wilson fermions explicitly. 

4.3.1 Naive fermions 

For r = 0, the meson integral equation becomes 



{m — AiaoGif 



(m — AiaoGif' 



$(x) 



X 



1 — X 



1 



(m — AiaQGif 



—7+ + 0:7 + (m - 4zcroG'i) + 




(4.22) 



2x(l - x) 



CHAPTER 4. THE MESON SPECTRUM 



45 



X 



dy 



(2vr) 



(m - AiaoGi] 
~2 2^ 



7+ - (1 - x)7 + (m - AiaoGi) 



27r 1 - X 



In this case, unbroken chiral symmetry of the action implies that the chiral limit of the 
theory is at m = 0. 



4.3.2 Wilson fermions 

For r = 1, the meson integral equation becomes 




+ 0:7 + m + 



9W 



271X 



2<!>{y) + J2rHy)Y 



27r(l -x) 



(4.23) 



This equation is quite close in structure to the 't Hooft equation, because the tadpole 
renormalization of quark mass vanishes. 

It is easy to see that the transverse part of the interaction kernel vanishes for the 
component $_;„j„2, for any value of k. It follows that ^--^712^^ an eigenstate of the equation 
with eigenvalue zero, when m = 0. This is the pseudoscalar Goldstone boson of the theory, 
with the corresponding wavefunction $i.„j„2 ~ J^^^ i'^ '^^^^ of the 't Hooft model. 
The fact that the chiral limit of the theory is at m = is a remarkable result, quite likely 
connected to the fact that the Wilson fermion transverse hopping term is constructed using 
projection operators in Dirac space. We feel that with a simple structure and an exact 
chiral limit at m = 0, numerical investigations of ()4.23|1 will not be too complicated. 

A result of particular interest is the Goldstone boson decay constant, since that allows 
us to estimate the overall scale of the theory. It is the same as that for the 't Hooft model, 
except for the change in dimensionality of Dirac space, 



{3+l)-dim 



V2(A)a+i 



1 



(l+l)-dim 



Using Nc = ?> and the experimental value ~ 
lattice cut-off, 7r/a_L ~ 300 MeV. 



(4.24) 

a_L V vr 

132 MeV, we obtain an estimate of the 



4.4 Outlook 



It is important to study bound states in QCD from first principles, because detailed quan- 
titative understanding of hadronic phenomena requires that. We have approached this 
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problem in the transverse lattice QCD formulation, exploiting properties of large-iVc and 
strong transverse gauge coupling limits. In these limits, we have explicitly obtained the 
generating functional of (3 + l)-dim QCD, starting from only quark and gluon degrees of 
freedom and no ad hoc model assumptions. Our study has been limited to mesons, and 
we have derived an integral equation satisfied by the meson light-front wavefunctions. Ex- 
traction of detailed properties (e.g. masses and form factors) from this equation requires 
numerical solutions to this integral equation. To study baryons, we have to find soliton 
solutions to the effective action calculated here in different baryon number sectors, and to 
study glueballs, we have to go to subleading orders in strong transverse gauge coupling 
expansion. 

Other groups have also studied bound states using the transverse lattice approach. We 
refer the reader to the detailed review of Burkardt and Dalley for a survey p]. They study 
transverse lattice QCD in the Hamiltonian formulation, using the colour-dielectric expan- 
sion. In the dielectric formulation, SU{N(.) group elements representing physical gluon 
fields, are replaced by general matrix variables corresponding to gluon fields smeared over 
short distance scales. These dielectric variables carry colour degrees of freedom, and form 
an effective gauge theory with a simple vacuum structure. The number of possible interac- 
tion terms in this effective gauge theory is much larger than in QCD. The requirements of 
symmetry restoration help in organization of these interaction terms, in a hierarchy suit- 
able for practical solutions. To date, all such solutions have been obtained in the large-Nc 
limit, which has restricted them to mesons and glueballs. Moreover, they have not been 
able to sum the transverse hopping expansion to all orders, and so they can not extract 
detailed chiral properties from their results. 



Appendix A 

Notations and Conventions 



A.l Space-time 



A point in (l+3)-dim Minkowski space time is specified by the coordinates x^, = 0, 1, 2, 3. 
x° is the time coordinate, while x^, , are spatial coordinates. The metric tensor of 
(1 + 3)-dim Minkowski space is 



/ +1 \ 
0-100 

0-10 

\ -1 / 



(A.l) 



and the invariant distance is 



ds'^ = g^^dx^dx" = {dx^y - {dx^y - {dx^y - (dx^) 



n2 



„2\2 



(A.2) 



It is implicitly understood that repeated indices are to be summed over their appropriate 
range, unless explicitly stated otherwise. The space-time indices /ijU, . . . are raised or 
lowered by contraction with the metric g'^'^ or g^^^,. 



V^ = {V',-V\-V^-V') 



(A.3) 



Derivatives with respect to contravariant (x^) or covariant (x^) coordinates are abbreviated 
as 



d., 



d 



The light-front coordinates are defined as 



_9_ 
dx„ 



: (X° ± X^) 
^ = 1= 

V2 



[Xq ± Xi) ± / 2 3\ I \ 

X± = -j^ , X X_L = (X2,X3), 



(A.4) 



(A.5) 
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where , x and are the light-front time, longitudinal and transverse coordinates 
respectively. The invariant distance in terms of light-front coordinates is 

ds^ = 2dx+dx- - {dx^Y - {dx^y . (A.6) 

For any four- vector V^, = while = —V±. In terms of the light-front components, 
the scalar product of two arbitrary four-vectors and takes the form 

V'Wi, = V+W- + V-W+ - V^W^ - V^W^ . (A.7) 

If is the four-momentum, then p~ , p+ and are the light-front energy, longitudinal 
momentum and transverse momentum respectively, = P/zP'' = for a real particle of 
mass m. The light cone energy p~ is then 

(p^Y + mF' 

P- = ^4^. (A.8) 

The (1 -|- 3)-dim transverse lattice space-time is the direct product of (1 -|- l)-dim 
Minkowski continuum and planar square lattice with lattice spacing a^. We denote the 
coordinates of a general point in this space-time by (x,x^), where x = x^ — {x^,x^) and 
x-^ = (x^.x'^) label the coordinates in (l-l-l)-dim Minkowski continuum and planar square 
lattice respectively. 

A. 2 Pauli matrices 

The Pauli sigma-matrices a^, and are the generators of the fundamental representa- 
tion of SU{2), i.e. the group of 2 x 2 unitary matrices, with unit determinant. 

They satisfy 

= ^^■'^ + ie*^■V^ (A. 10) 

The raising and lowering operators are c"^ = | (c^ ± icr^), 

.-(° (A.n) 
A. 3 Dirac matrices 

The Dirac gamma-matrices {7''}, with /i e {0, 1, 2, 3}, form a 4-dimensional representation 
of the Clifford algebra defined by the anticommumtation relations 

[7^71+ = 25^^^ (A.12) 
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Further define 

7' = 75 = ^7V7V- (A.13) 
From the anticommutation relations, it follows that 

[7'>7l+ = 0, {rf = l. (A.14) 
We use the particular matrix representation provided by 

In this case, 7^,7^ are hermitian, while 7^,7^,7^ are anti-hermitian. Also 

^0(^M)t^O ^^M. (A. 16) 

On the light-front, 

7± = -^(7°±7^)- (A.17) 
It follows that (7''")^ = (7")^ — 0) while 7"^7~ + 7~7"'" = 2. Explicitly, 

v/2cT+ \ _ / V2(r- 



A.4 Lie algebra of SU (Nc) 

The Lie group SU (Nc) is the group of A^c x matrices of unit determinant. Every element 
V of SU {Nc) can be parametrized as in terms of — 1 real parameters labeled as cua 

V ^expliY^taCUaj , (A.19) 

where ta are the generators of the Lie algebra of SU (Ac) and cUa are real parameters. The 
generators are traceless, Hermitean Ac x A^c matrices, i.e. 

tr {ta) = , tl = ta. (A.20) 

They satisfy the commutation relations, 

[ta,tb]_ = ifabctc ■ (A.21) 
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where the structure constants fabc are completely antisymmetric and real. They also obey 
a completeness relation 

a=l ^ 

We use the normalization convention 

tr (tah) = Sab ■ (A.23) 

For a general matrix M in both colour and Dirac space, we denote its elements as M^^ 
(Latin indices for colour and Greek indices for spinor). For its trace, we use the notation 

• trace in colour space: tr(M), 

• trace in Dirac space: trD(M), 



trace in colour and Dirac space: Tr(M). 



A. 5 Parallel transporters 

In a gauge theory, parallel transporters are unitary matrices, U (Cy-x), forming a map from 
directed curves Cy-x to the gauge group. They satisfy the group properties: 

• U (0) — 1 where is a curve of zero length. 

• U (C2 o Ci) = U (C2) U (Ci), where C2 o Ci is the path composed of Ci followed by 
C2. 

• U (— C) = U {C)~^ = U {C)\ where — C denotes the path C traversed in the opposite 
direction. 

Under a gauge transformation specified by V{z), 

U {Cy.,x) ^ C/' {Cy.,x) = V{y)U {Cy.,x) V'' (x) . (A.24) 
Given the gauge field (x), for an infinitesimal curve, 

U {Cx+dx;x) ^1-A^{x)dx'' + ---. (A.25) 
If c'^{t) parametrizes the curve C from c'^(O) = x'^ to c'^(s) = y'^, then 



d fir'' 
With the initial condition (0) = 1, the solution is 



U (C) (A.26) 

t=s 



C/(C)=7'exp|-^ A^{c{t))^dt^ =Veicpi^- j A'' {x,x^) dx^^ (A.27) 
where the symbol V denotes path ordering with respect to the variable t. 



Appendix B 

The 't Hooft Model 



B.l Preliminaries 

The 't Hooft model is the (l + l)-dim SU (Nc) gauge theory, with quarks in the fundamental 
representation, in the large-iVc limit. We describe various aspects of its solution in this 
Appendix. We use functional integral methods, since that allows the large- Ac limit to be 
imposed in a straightforward fashion. 

The action for this theory in terms of the quark and the gluon fields is 

Nf 

S=^J d'xtr{F>^-'{x)F,,{x))+ J d'xjya{x)[trD,-ma]^a{x). (B.l) 

When the bare quark masses rUa are zero, then the action is invariant under: 

• Single flavour case: Uv{l) ^7a(1)- 

• Multi-flavour case: -5C/l(A/) ® SUR{Nf) ® C/y(l) ® Ua{1). 



B.2 Generating functional for quark- ant iquark 
correlation functions 

In presence of an external quark bilinear source term, the generating functional of the 
theory is (for sake of simplicity, we work with a single flavour), 

Z[J]= J [D^ • DW\ [DA] exp [il A] + (gauge terms)] , (B.2) 

I-S + J d'xd'y^.{x)J^^{x,y)^^{y) . (B.3) 
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Here "gauge terms" stands for both the chosen gauge as well as the corresponding ghost 
determinant. Arbitrary connected time-ordered quark-antiquark correlation functions can 
be obtained by taking functional derivatives of W [J] = —ilnZ [J] with respect to J, e.g. 



T < (x) Vl/f (y) 



6W [J] 



(B.4) 



j=o 



6^W [J] 



5Jif {x,y)5Jli {z,w) 



(B.5) 



j=o 



In order to compute the quark propagator and the quark-antiquark scattering amplitude, 
W[J] should be known to 0(7^). 



B.3 Computing Z [J] in large- A^c limit 



We work in the ghost-free light-front gauge = A'^ + A^ = 0. In this gauge, the functional 
integral over the A~ variables is purely Gaussian; so A~ can be integrated out exactly. 



Z[J] 



[D-^ ■ DW\ [da-] exp 
X exp 
X exp 



d'xd'y ^:{x)J^^ix,y)^^{y) 



(B.6) 



d X -^ij (^) 



-6u6jkd'iA,^i [x] 



d'x (<(x)(7+)"'^f(a:))Ar.(x) 



The integration over gauge fields produces a non-local interaction term which is quadratic 
in fermion bilinears. Suppressing the overall normalization constant. 



Z[J] 



[D^ ■ DW\ exp 



t / rfWi/<(x)J"^(x,y)v&f(y) 



X exp 
X exp 



I / d^x ^"{x)6,, (i {Yf^d^ - "^5"^) ^i(^) 



(B.7) 



d'xd^ h{x-y) (7+)"' (7+)"'' <(a;)^f (y)^(i/)vi/5(x) 



7/3 -j-a. 



where h{x — y) is the solution of d^h{x — y) = 5^'^\x — y), with the principal value 
prescription as the boundary condition. To carry out the fermionic integration, we use the 
identity 
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I_[UadCadQexpTTf{{^,U}) 
= I [Ua diadla] [llab ddab^ ((7a6 - la^h)] exp Tl/ {{(J^d}) 

= / [ria diadl^ [Hab daabdXab] exp [iAafe (daft - C^)] exp Tr/ ({cTcd}) (B.8) 

= / Wlab daabdXab] dct exp [iAafeCTab + Tl/ {{(Jed))] 

= I lHab dcTabdXab] exp (Ti [ill (-^A) + i(T^A + / ({o-cd})] ) 

This identity allows the functional integral to be written only in terms of bilocal bosonic 
variables: 

(B.9) 



Z[J]= J [D'J ■ DX] exp \iV {a, A; J)] , 



where 



Via,X; J) 



d^x Tr [ln(-iA(a;, x))]+ / d^xd^y X'^f {x, y) aff {x, y) 



a.3 



+ (x, y) [5,, {z {^^r^ d, - m5«^) 5^^) {x - y) + J^f (x, y) } (B.IO) 



h{x-y) (7+)"' {^^y'af{x,y)a]Uy.x) 



2Nr 



V (cr, A; J) is of O (iVc) , so the large-iVc limit of the functional integral amounts to evaluating 
the functional integral at its stationary point. (This illustrates why the Nc ^ 00 limit can 
be thought of as a classical limit of a quantum theory; albeit different than the classical 
limit corresponding to — > 0.) The equations determining the stationary point are: 



dV 



^Ki i.^^ y) 

dV 



dcTif {x, y) 



<j,X 



<T,\ 



i[x-')yy,x) + arf{x,y)^0 (B.ll) 

Xl^ {x, y) + 5,, (i {^^r^ d, - m5«^) 5^'^ {x-y) + (x, y) 

(B.12) 



-|r^. M" {l^rh{x-y)ali{y^x)=^ 



The first stationary point equation implies 

Xj{x,y) =i{a-^f." {y,x) . 
Substitution of this in to the second stationary point equation yields 

i (a-^) J; (y, x) + % {i irr' d, - m5"^) 5(2) {x-y) + jff (x, y) 



(B.13) 



,7/3 



— 75 



(B.14) 
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Solution of (1B.14|1 will yield a as a function of J. Thus we have 

Z [J] = exp \tV {a, A; J) ] W [J] = V (a, A; j) = V,s (a; J) , 

where Ves (o"; J) , obtained by substituting for A in terms of a, is given by 



(B.15) 



Kff (cr; J) ~ '^^3; Tr [In cr(x, a;)] + J d^xdPy 

'<! V) (x, y) + 5., (i {^r^ d, - mS^^) 5^'^ (x-y)] (B.16) 



h{x - y) {-i^Y^ {l^y^ erf (x, y) a]j {y, x) 



2Nr 



— 75 



Note that V^s (o"; J) depends on J explicitly, as well as implicitly through a. 

The generator of IPI vertex functions T [{p] is the Legendre transform of In 
terms of the effective field ip conjugate to the external source J, 



sw [J] 



(B.17) 



W[J]- j d'xd'yvf{x,y)Jtl'{x,y) 
i J (Px Tt [\nLp(x, x)] 

+ j d'xd'y [iftf {x, y) 5,, (z (7'^)"'' d, - m(5°^) 5^'^ {x - y) 



(B.18) 



h{x - y) (7+)"^ (7^)^^^ '^f {x, y) {y, x) 



7(5 



2N, 



It can be easily seen that in the large- Ai'c limit, (p = a + O (l/Nc). Also 



[1^(2)] 



a/3;7<5 
ij;kl 



[x,y; z,w) 



ij;kl 



[x,y;z,w) 



(B.19) 



where W^'^^ and T^"^^ are functional second derivatives of W[J] and T[lp] respectively. 



B.4 Quark propagator 



Let us first solve the stationary point equation ()B.14p when J = 0. 

- (ao^) (y, x) + i5,, (i {^T^ d, - m^^) 5^'^ {x - y) 



■ \ 75 



N. 



(B.20) 
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(Jo is expected to be invariant under space-time translations, and so can be parametrized 
as 



(CTO 



x,y) 



(2n 



1 



m — S (p) + ze 



/3a 



-ip-(y-x) 



(B.21) 



With this parametrization, the J = stationary point equation becomes an integral equa- 
tion for S (p) , 



(2^ 



1 



7 



1 



^ — m — E (g) -|- ie 



7^ 



_ q+y 

where the principal value prescription regularizing the singularity is defined as 

1 



(B.22) 



' 1 " 


1 


1 




~ 2 





.z - le) 



(B.23) 



(IB.22|1 is just the 't Hooft integral equation for the quark self-energy. The Dirac structure 
of the equation requires E (p) to be proportional to 7"*". The instantaneous nature of the 
kernel demands that E {p) depend only on p"*". The equation can then be solved for E (p), 



E(p) 



27rp^ 



(B.24) 



From (lR4ll . (1051) . ^^\^ and (jR2nll . it follows that the J = stationary point, ao, is 
the quark propagator. Explicitly, 



—I 



—I 



{2tx 
(2 



1 0a 



I <7^7^ I 



27rp+ 



e-ip-(y-^) (B.25) 



p 



■2 — + — + 



/3a 



-ip-{y-x) 



B.5 Meson spectrum 



The time-ordered four-point quark-antiquark Green's function can be obtained from either 
W [J] or r [if]. Combining (Ib311 and (IHToll . 

T(^r(x)vl/^^(2/)^(.)vl/fH)) - (i:{W:{x)^>'^{y))){T{%{z)¥,{w))) 

6^W [J] 



^Jtf {x,y)SJ2i {z,w) 



j=o 



(r(2))-^ 



a/3;75 
ij:kl 



{x,y;z,w} 



(B.26) 



1/3=0-0 
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Using r [ip] in (IB.18|) . one can show that 

r(2)| 



I l/3=CT0 



where 



S^M^ {x, y, z, w) = -i ((To {y, z) [a^ ^)"" [w, x) 



S + K 



l\/37 



(B.27) 



(B.28) 



K^^' {x, y; z, w) = (7"^) (7"')''^ Hx - y)6^'^ {z - y) 6^^ {w - x) . (B.29) 

If there are mesons in the spectrum, they will show up as poles in the four-point quark- 
antiquark Green's function. To find the positions of these poles, all we need to do is find 
the zeros of [5* + K]. The meson spectrum is thus determined by 



(B.30) 



Transforming this equation to momentum space, one gets the integral equation 



— m + 



9 4- \ Ct<5 

9 1 



52 - „ + 



27rj9^ 



d kid k2 



■ p 



xt ih, h) = o. 



(B.31) 



{pt-kty 

Without loss of generality, we scale the momenta and choose 

Pi = {pt^Pi) > P2 = (P^,P2 ) , P| = 1, p| = X, kf = y = yi, k^ = y2. (B.32) 
Then the momentum space integral equation ()B.3H1 becomes 

»i + + m + ^ 

xtf (1 -a;,P2;a;,pr) 



9 1^ 
2nx 



X 



X — + ^ + ie 
-P2 7+ - (1 - x) 7- + m - 



<5o 



2p2 (1 ~ 3;) — + ^ + 



(B.33) 



X 5(1 -7/1 -1/2) P 



,5 f dk, dk^ dyidy^ . . 

2 <J \Pi ~r P2 ~ 1^1 ~ 1^2 ) 



The light-front wavefunction is defined through the projection 



x7 (1/1,^2 ) 



(x) = y dp-^ dp2 6 - - P2 ^ x"/ (1 - x,^j2 ; ) . (B.34) 
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The content of (|B.33j) can then be written solely in terms of (x) as, 



m 



m 



X 



1 — X 



(x7 + m)'^^ (— (1 — a;) 7 + m) 
2x{l-x) 



[X 



5a 



(B.35) 



^'^^ {y) ■ 



Using standard matrix notation, ()B.35|1 can be written in a compact form, 



M^^ix) 



X 



+ 



1 — a; 



$ ix) 



(B.36) 



+ 



A-nxil — x) 



dyP 



[x'y + m) 7"'"$(?/)7^ (— (1 — x)7 + m) 



Until now we were working with a single flavour. The calculation can easily be extended 
to the multi-flavour case, since quark-antiquark pair production/annihilation vanishes in 
the large- A^'c limit. If and mb are the masses of flavour a and h respectively, then the 
multi-flavour analogue of ()B.36|) is 



M'^^{x) 



ml -31 



x 



+ 



TT 



1 — X 



(B.37) 



+ 



Attx{1 — x) 



dyP 



_{x - yf 



[x-1 + ma) 7+$(y)7+ (-(1 - a;)7 + rrn) 



From the structure of above equation, one can see that only one component of ^"'^ (x) is 
independent. We expand $ as 



$ (x) = $1 (x) 1 + $+ (x) 7+ + <l>^ (x) 7 + $+_ (x) - [7+, 7 ] _ • 



(B.38) 



Substituting this expansion in (IB.37jl . one can easily see that $_ (x) satisfies the 't Hooft 
equation, 

r,2 ' 



H^ix) = M2$„ (x) 



ml -31 



X 



mt - 



1 — X 



$_ (x) 



TT 



dyP 



(y) 



(B.39) 



The remaining components $1 (x), $+ (x) and $+_ (x) are algebraically determined in 
terms of $_ fx) as 



^ , . I f rua mt, 
/ V X 1 — X 



$_ (x) 



(B.40) 
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(X) = + (X) . (B.42) 

B.6 't Hooft solution 

Over the years, the properties and solutions of the 't Hooft equation ()B.39j) have been 
investigated in great detail. Prominent results are [26] [TO] [T2] : 



For finite norm solutions, must vanish at the boundaries x = and x = 1 as 

x^'^ and (1 — x)^*" respectively, where 

2 2 

rrPa cot (nPa) = 1 - ^ , nP, cot (nP,) = 1 - !^ . (B.43) 

The spectrum of is purely discrete. The light-front wavefunctions $"(x) are 
gauge invariant, and they form a complete and orthogonal set. 

^- i^W- i^') = - x') (B.44) 



rfx<l>!!*(x)<l>'!(x) =5"'^ (B.45) 







When rria = nib = 0, the lowest meson state has zero energy. The corresponding 
wavefunction is $i(x) = 1. This is the Goldstone boson corresponding to chiral 
symmetry breaking. 

For large n (large energy), the eigenvalues and eigenfunctions can be approximated 
as 

-Tsn/N /;^ / ,U>1. B.46 
$" (x) ~ V2sm (vmx) J ^ ^ 

Thus the eigenvalues form an almost linear trajectory in the — n plane. 
Using the operator K defined by 

ir<l>_(x) = P / dy—^, (B.47) 

Jo [y ~ ^) 

and the fact that ^ 

[ dx (x) [H, K] $'! (x) = , (B.48) 
Jo 
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one can prove that 



/ ay- 



1 $-(^) 

y 



(B.49) 



Here P„ = (—1)" is the parity of the eigenstate n. The lightest meson is a pseu- 
doscalar, Pi = —1. 

Continuous symmetries cannot be spontaneously broken in two dimensions. In the 
't Hooft model, chiral symmetry is realized in the Berezinskii-Kosterlitz-Thouless 
mode EH, 



*^(0)) ~ x"^/^^ . (B.50) 

To obtain results for the broken chiral symmetry phase (e.g. non-zero chiral conden- 
sate), one must take the ^ oo limit before taking the m — limit. 

The chiral condensate, regularized by a split-point definition and by subtracting the 
free field theory contribution, is [7] 



-iVc 



m 



27r 
1 - 



1 + 7i5 + In 



m 



9' 



n=2 



1 1 

Pn n 



(B.51) 



where /5„ G [ra— 1, are the solutions to the 't Hooft boundary condition 7r/3 cot (vr/?) 
1 — (m'^n/g^). In the chiral limit, 



9' 



12 V TT 



(B.52) 



Mesonic matrix elements of quark-antiquark bilinear operators are easily obtained 
from ()B.4nilB.42|l [IQJ. In particular, the axial current coupling is saturated by the 
Goldstone boson in the chiral limit, and the pion decay constant is given by 



(B.53) 



When one of the quark flavours is much heavier than the other, Mq ^ m, one can 
take the static limit of the 't Hooft equation 0. In terms of the non-relativistic 
variables, E^tat = M - niQ, t = {1 - x)mQ and X-{t) 



^$_(1 - ^1 



£^statX-(^) 



2t 



t 

+ 2 



dyP 

Jo 



X-is) ■ (B.54) 



Detailed results for physical properties have to be obtained numerically. 



Appendix C 



Lattice QCD at strong coupling and 
large- A^c 

C.l Preliminaries 

The d-dim Euclidean lattice QCD action in the strong coupling limit is 

Se = - ^ \^{x + ij){r - i-f^") U (x, /i)^(x) + '^{x) (r + ^7^) U\x, + /i) 

+im^'*(a;)^(x) . (C.l) 

Here x = (x^x^, . . . ,x'^) are points on a d-dimensional hypercubic lattice, and the lattice 
directions are labeled hj fi = 1,2, . . . d. We choose to work with even d, and use units such 
that the lattice spacing a = 1. The Dirac matrices have dimensions C = 2^^/^, and satisfy 

[7^71+ = -2^^'^ , {in^ = -r. (C.2) 

The gauge links U{x, fi) connect the nearest neighbour points on the lattice. They are 
independent random variables in the strong coupling limit (there is no pure gauge term 
in the action). The fermion nearest neighbour hopping term is written using the Wilson 
parameter r. With naive fermions (r = 0), the theory gives rise to 2*^ fermion zero modes, 
one for each corner of the Brillouin zone. Choosing r 7^ removes all fermion doublers, 
leaving behind only the fermion mode corresponding to the continuum theory. 
For naive fermions, (IC.ljl reduces to 

Sen = -^^['^(x + /i)7''f/(x,/i)^(x) -'*(x)7^?7"^(x,/i)^(x + /i) 

X,fJ. 

+im^*'(x)^(x) (C.3) 
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The global symmetry group of this naive fermion action, Sen, can be easily deduced by 
rewriting it in terms of spin-diagonalized fermion variables x ^-nd x- Let 

= A{x)x{x) , = x(^)At(x) . (C.4) 

The choice ^ 

A{x) = {j'f {^Y^ ■ ■ . {YY , A\x)Aix) = 1 , (C.5) 

reduces the 7-matrices in the action to position dependent phase factors, 



A^{x + fi)j''A{x) = A^'{x) ■ 1 



(C.6) 



where 

A'^(x) = (-1)^''<M^\ (C.7) 

Thus the action reduces to 



'Ejv = ~z2^ — \x{x + i^)U{x,ii)x{x) +x{x)U^{x,ii)x{x + 11) 

+im^x{x)x{x) . (C.8) 



It is easy to see from this form of the action that the fermion hopping term is invariant 
under Uo{C) ® Ue{C) ("o" and "e" respectively stand for odd and even points): 



X (X) = Uo(e) X{X) \ r A ..^ X 

_,, . -/ wrt r for > odd (even) (C.9) 

X'{x) = x{x) t/e(o) J \^ J 



X{x) X'{x) = f/o(e) X{x) 



This invariance group describes the chiral symmetry of naive lattice fermions, and it is 
considerably larger than the corresponding chiral invariance group in the d-dim continuum 
theory. The fermion mass term is invariant only the diagonal subgroup Ua{C) of Uo{C) <S> 
Ue{C), defined by Uo = Ue. 

For fermions with r 7^ 0, the action Se is invariant only under Ua{C), even if m = 0. 
The r-dependent terms break the non-diagonal part of Uo{C) ® Ue{C) explicitly. 

C.2 Generating functional for quark-antiquark 
correlation functions 

Since we are interested in meson correlation functions, we add an appropriate meson source 
term to the starting action. The partition function of the theory is 



Z [J] = J [D^ ■ [DU] exp 



(C.IO) 
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The integral over the gauge links U{x, n) completely factorizes to a separate integral for 
each individual link, and so can be evaluated explicitly in the large- A^c limit. The result is 



Z[J] 



[D^^ ■ DW] exp 



X exp 



X,fl 



where 



f[B] 



1 + —B - 1 - In 



11 1 



Matrix elements of B{x,fi) in Dirac space are 



(C.ll) 



(C.12) 



(C.13) 



The integration over gauge link has produced non-local interactions between fermion bi- 
linears. To perform the fermionic functional integral in the large limit, we first rewrite 
the fermionic integral in terms of bosonic variables using the identity ()B.8|) : 



Z[J] = j [Da ■ D\] exp [-V{a, A)] 



(C.14) 



V{a,\) = - J]Tr[ln(-zA® l,)(a;)] - J]2(t"^(x)A"^(x) 

X X 



Here 



/[a] = 5^Tr/[fi(x,/x)], 

X,fl 

with matrix elements of B°'^{x,^) in terms of a given by 

B-^x^ ^^) = -a^ix + /.) (r + Z7j)'^ (r - ^7?)^' 



(C.15) 



(C.16) 



(C.17) 



In the large- Ai'c limit, the bosonized functional integral reduces to its stationary point value. 
The stationary point equations are: 



dV 



_ = -ta^P{x) - N, (^A (x) = , 

7, A 



(C.18) 
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dV 



-Nc 5^ [ (r + nj) ^{x - /i) (r - Z7j) /' [B^x - /i, ^)] 



Here 



f[B] 



/3a 



Solution of the first stationary point equation ()C.17jl gives 



l\Pa 



. (C.19) 
(C.20) 

(C.21) 



Substituting this in the second stationary point equation ()C.18jl yields a self-consistent 
equation for a, 

Nc {a-^f" (x) + {im6''^ - J"^(x)) 



+ ^7j) /' [B{x, /x)] a{x + fi) {r + 27^) 



13a 



. 



(C.22) 



Solution of this equation will yield cr as a function of J. The large-iVc generating functional 
of connected meson correlations is thus given by 



W[J] = -\nZ [J] = V {a, A) = Ks (a; J) , 
with Ves (o"; J) (obtained by substituting A in terms of a ) given by, 

Vesi'^;J) = ^Tr[\n{a^U)ix)]+im^tiY)a{x) -^tij:)j'^{x)W{ 



(C.23) 



(C.24) 



X,fJ, 



B{x,fi) 



2 2 



The matrix B is B, eq. ()C.17|) . evaluated at the stationary point, 



(C.25) 



The generator of IPI vertex functions T [ip] is the Legendre transform of [J]. In 
terms of the effective field ip conjugate to the external source J, 



cp'^^ix) 



dW [J] 
'dJ^Hx) 



(C.26) 
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r 



with 



X 

Tr [In {if ® Ic)] (x) + im ^ trDV5( 

X 

1 



X,fJ, 



iV2 



Note that 9? ^ cr in the large- A^c limit. Also 



(C.27) 



(C.28) 



(C.29) 



where W'-^^ and T^^) are functional second derivatives of W[J] and r[ip] respectively. 



C.3 Chiral condensate 

When J = 0, we expect the stationary point to be translationally invariant. Then equation 

rf(l-r2)ao 



221) gives, 



— h zm H 

O-Q 1 



/I + (1 _ r^) ^2 







where o"o^(a;) = Ai'ccro^"^. This equation can be solved to yield 



(Jo = t 



-m{d - 1) + rf^m2 + (2d - 1) (1 - r^) 



(i^ (1 — r^) + 



The chiral condensate is given by 

In particular, 

• Naive fermions (r = 0): 



dW [J] 



a J"" fx) 



(ao)"" (x) = N^Cao 



(C.30) 



(C.31) 



(C.32) 



j=o 



-m 



{d-l) + dVm^ + 2d-l 



(C.33) 



In the chiral limit m — 0, a non-zero chiral condensate signals spontaneous break- 
down of UoiC) ® UciC) to its diagonal subgroup f/A(C), 



lim {mm) = iN^C 



V2d-l 



d 



Wilson fermions (r = 1) 



m 



(C.34) 
(C.35) 



Chiral symmetry is explicitly broken in this case. 
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C.4 Meson propagator 

The meson propagator is given by 



<(a;)^f(a;)^(l/)v^5(l/) 



dW [J] 



dW [J] 



+ 



d^W [J] 



j=o 



dJ''i^{x)dJ^^{y) 



j=o 
(C.36) 



(/3=(T0 



where 



ip=cro 



{G, + (1 - r') a',G,) 5^ { (r - ^J,r' (r + ^J,y^ 6^.,, 



+ {r + i-ff,)"' (r - i-fi.y'^ 



with 



-2d (1 -r^) a^Gsrt''^^. 
1 



x,y 



2 1 + v/l + (l-r2) 



^2 



-1 



4^1 + (l_r2)a2 (l + ^l + (l_r2)a2 



(C.37) 
(C.38) 
(C.39) 



Eq. (l(y.37|l is the inverse of the connected meson propagator. The poles of the meson 
propagator therefore correspond to the zero eigenvalues of r'^^)(p). Transforming T^'^\x,y) 
to momentum space, we get 



(^-2rf(l-r2)'a^G2)rt^^ (C.40) 
+ (Gi + (1 - r') a^G,) J2 {2 + 7,"S?) cos (p,) 

+2r {6"'jf - 7^5^^) sin (p,) 



T^'^\p) is a matrix in the G^-dimensional space of the Clifford algebra generated by the 
d gamma matrices. This matrix has to be diagonalized in order to find the spectrum of 
meson states. In 4-dim, the commonly chosen Clifford algebra basis is: 



Tm = {5 = 1, 1/ = 7^, T = -[7^, 7^], A = 7^75, P = 75} . 



(C.41) 



There is at most one meson state for each spin-parity sector, i.e. no radial excitations 
survive in the strong coupling limit. 
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For naive fermions (r = 0), the inverse connected meson propagator is 

(r(2))"^'-^ (p) = -NMl^-2dalG,)5-^5-^+2 [G, + alG,) ^.'if cos (p,) ] , (C.42) 

L \(Tq / J 

and it is exactly diagonalized by the basis ()C.4H) . The resultant dispersion relation is 



where n is (pseudoscalar), 1 (vector), ... , c? — 1 (axial), d (scalar), and is +1(— 1) 
when Tm anticommutes (commutes) with 7^. This result shows that, the fermion doubling 
phenomenon for naive fermions, produces C^ pseudo-Goldstone bosons with mass 



= -77.^^^ = .m^^^ , (C.44) 
ao (Gi + (T^Ga) (To 

one for each corner of the Brillouin zone. The pseudoscalar located at p = is identified 
as the physical pion, while the other modes at p = (i.e. n 7^ 0) are interpreted as the 
heavier meson states [TH]|19|. 

For Wilson fermions, we restrict ourselves to li = 4. For r 7^ 0, the inverse connected 
meson propagator separates into three block matrices corresponding to S, PA and VT 
channels. As r increases from to 1, states within each channel get more and more 
strongly coupled. For r 7^ 0, the "chiral limit" is non-perturbatively defined as the value 
m = rric for which the pseudoscalar mass vanishes: 

r — > : TTic ~ ^ r , r = 1 : iric = 2 . (C.45) 

At m = mc, (To decreases smoothly from ^ to | as r increases from to 1. For r 7^ 0, 
all physical meson states are located at p = 0. At r = 1 in particular, there are only four 
meson states, one pseudoscalar and three vectors [nj[15j[20]. 
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